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Abstract 

We propose a model for (unilateral) contact with adhesion between a viscoelastic body and a rigid 
support, encompassing thermal and frictional effects. Following Fremond's approach, adhesion 
is described in terms of a surface damage parameter X . The related equations are the momentum 
balance for the vector of small displacements, and parabolic-type evolution equations for X and 
for the absolute temperatures of the body and of the adhesive substance on the contact surface. 
All of the constraints on the internal variables, as well as the contact and the friction conditions, 
are rendered by means of subdifferential operators. Furthermore, the temperature equations, 
derived from an entropy balance law, feature singular functions. Therefore, the resulting PDE 
system has a highly nonlinear character. 

The main result of the paper states the existence of global-in-time solutions to the associated 
Cauchy problem. It is proved by passing to the limit in a carefully tailored approximate problem, 
via variational techniques. 

Key words: contact, adhesion, friction, thcrmoviscoclasticity, entropy balance. 
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1 Introduction 

In this paper we propose and analyze a PDE system modelling thermal effects in adhesive contact 
with friction. More specifically, we focus on the phenomenon of contact between a thermoviscoelastic 
body and a rigid support. We suppose that the body adheres to the support on a part of its boundary. 

Following M. Fremond's approach [551 HH]: we describe adhesion in terms of a surface 
damage parameter X (and its gradient VX), related to the state of the bonds responsible for the 
adherence of the body to the support. Further, we consider small displacements u and possibly 
different temperatures in the body and on the contact surface. Internal constraints, such as unilateral 
conditions, are ensured by the presence of non-smooth monotone operators, also generalizing the 
Signorini conditions for unilateral contact to the case when adhesion is active. Frictional effects are 
encompassed in the model through a regularization of the well-known Coulomb law, here generalized 
to the case of adhesive contact and assuming thermal dependence of the friction coefficient. 

The mathematics of (unilateral) contact problems, possibly with friction, has received notable 
attention lately, as attested, among others, by the recent monographs on the topic [22] and [3T]. The 
analysis of Fremond's model for contact with adhesion has been developed over the last years 
in a series of papers: in [7]-[H] we have focused on existence, uniqueness, and long-time behavior 
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results on the isothermal system for (frictionless) adhesive contact. Thermal effects both on the 
contact surface, and in the bulk domain, have been included in the later papers [9]-|10j. respectively 
tackling the well-posedness and large-time behavior analysis. In we have dealt with a (isothermal) 
unilateral contact problem taking into account both friction and adhesion. As far as we know, the 
present contribution is the first one addressing a model which combines adhesion, frictional contact, 
and thermal effects. 

The PDE system. Before stating the PDE system under investigation, let us fix the notation for 
the normal and tangential component of vectors and tensors that will be used in what follows. 

Notation 1.1 Given a vector v G R 3 , we denote by «n and vt its normal component and its 
tangential part, defined on the contact surface by 

un := v • n, vt := v - %n, (1.1) 

where n denotes the outward unit normal vector to the boundary. Analogously, the normal com- 
ponent and the tangential part of the stress tensor cr arc denoted by on and cr t , and defined 
by 

ctn := crn • n, er T := crn <r^n. (1-2) 

We address the PDE system 
d t (ln(tf))-div(d 4 u)-Atf = /i infix (0,T), (1.3) 

<W=<[° in(3Q\r c )x ( 0,n 

\ -k(X)(d - 0.) - c'(tf - 0.)W(-oo,o] K)|5 t u T | in r c x (0, T), 
9 t (ln(^))-5 t (A(X))-A^ = fe(X)(#-0 s ) + c'^-^)5/ ( _ OOi o]("N)|5tU T | in T c x (0,T), (1.5) 

d n J s = in dT c x (0,T), (1.6) 

- diver = f with cr = K e e(u) + K v e(d t u) + til in Q x (0, T), (1.7) 

u = inrix(0,T), crn = g inr 2 x(0,T), (1.8) 

cr N G -Xu N - 97(_ OOi0 ] (« N ) in T c x (0, T), (1.9) 
<t t G -Xu T -c(tf-tf s )ar ( _^ i0] (u N )d(d t u) inr c x(0,T), (1.10) 

9 t X + Sdl^M 9 **) ~ AX + M[o,i]W + ^PO 3 ~A'(X)(i? a - tf cq ) - ||u| 2 in T c x (0,T), (1.11) 
9 ns X = in<9r c x(0,T) (1.12) 

with 8 > , where 

if v T ^ 

r 1.13 

6 5i} ifv T = 0, 
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and B\ is the closed unit ball in Mr . Note that d : Br =4 R is the subdifferential of the function 
j : R 3 — >• [0,+cxd) defined by j(v) = |vt|, cf. (|2.5p later on. In what follows, f2 is a (sufficiently 
smooth) bounded domain in R 3 , in which the body is located, with dil = T± U T2 U T c and T c 
the contact surface between the body and the rigid support. From now on, we will suppose that 
r c is a smooth bounded domain of R 2 (one may think of a flat surface) , and we denote by n s the 
outward unit normal vector to dF c . Let us point out that u is the vector of small displacements, cr 
the stress tensor, whereas X is the so-called adhesion parameter, which denotes the fraction of active 
microscopic bonds on the contact surface. As for the thermal variables, 1} is the absolute temperature 
of the body fi whereas i? s is the absolute temperature of the adhesive substance on T c . Here and 
in what follows, we shall write v, in place of v\r c , for the trace on T c of a function v defined in fi. 

While postponing to Section [2] the rigorous derivation of the PDE system (|1. 3111. 12"]) based on 
the laws of Thermomcchanics, let us briefly comment on the nonlinearities therein involved. First, we 
observe that the singular terms ha(#) and ln(i9 s ) in (|1.3|) and (|1.5j) force the temperatures 1} and -d s 



2 



to remain strictly positive, which is necessary to get thermodynamical consistency. The multivalued 
operator die ■ K =4 K (with C the interval [0, 1] or the half-line (— oo, 0]) ) is the subdifferential (in 
the sense of convex analysis) of the indicator function of the convex set C . We recall that die (y) ^ 
if and only if y £ C, with dlc(y) = {0} if y is in the interior of C, while dlc(y) is the normal 
cone to the boundary of C if y € dC . By means of indicator functions, we enforce the unilateral 
condition % < , as well as the constraint X 6 [0, 1] and, if S > , the irreversibility of the damage 
process, viz. dfX < 0. Then, the subdifferential terms occurring in (|1.4p , (|1.5p . (|1.9p . (|1.10[) . and 
(jl.lip represent internal forces which activate to prevent «n , X, and dtX from taking values outside 
the physically admissible range. Finally, K e and K v are positive-definite tensors, k, A, ct, and c 
are sufficiently smooth functions, whereas "& oq is a critical temperature and h, f , and g are given 
data: we refer to Section 2 for all details on their physical meaning, and to Section I3T21 for the precise 
assumptions on them. 

As for conditions (|1.9|) - (|1.10|) . let us observe that (|1.9p can be rephrased as 



which, in the case X = , reduce to the classical Signorini conditions for unilateral contact. Con- 
versely, when < X < 1 , ctn can be positive, namely the action of the adhesive substance on T c 
prevents separation when a tension is applied. Moreover, in view of (jl.9p . (|1 . 10[) can be expressed by 



|cr T + Xut| = c(i? — #s)|cn + Xun| =>■ 3 v > : 9 t ux = — v(ai + Xut), 

which generalize the dry friction Coulomb law, to the case when adhesion effects are taken into 
account. Note that the positive function c in (|1.15p is the friction coefficient. 

Related literature. We refer to [32] and [3T], for a general survey of the literature on models 
of (unilateral) contact, as well as to the references in [8] for isothermal models of contact with 
adhesion, possibly with friction (see also [H]). Here we will just focus on (a partial review of) 
temperature-dependent models for frictional contact: as previously said, to the best of our knowledge 
the PDE system (|1 .311 1 . 1 2|) is the first one modelling adhesive contact with frictional and thermal 
effects. 

The major difficulty in the analysis of unilateral contact problems with friction is the presence 
of constraints on both the (normal) displacement «n and the tangential velocity <9tUx . It can be 
overcome only by resorting to suitable simplifications in the related PDE systems. Over the years, 
several options have bee explored in this direction, without affecting the physical consistency of the 
underlying models. 

Since Duvaut's pioneering work |20j . a commonly accepted approximation of the dry friction 
Coulomb law (|1.15p . which we will also adopt, involves the usage of a nonlocal regularizing operator 
31 , cf . (|1.16p below. An alternative rcgularization of the classical Coulomb law is the so-called SJK- 
Coulomb law of friction, which is for example considered in |28j and in pQ, dealing with quasistatic 
models for thermoviscoelastic contact with friction. 

An alternative possibility is to replace the Signorini contact conditions Q1.14p with a normal 
compliance condition, which allows for the intcrpenctration of the surface asperities and thus for 
dispensing with the unilateral constraint on mn ■ Analytically, the normal compliance law corresponds 
to a penalization of the subdifferential operator 9/(_oo.o] in (|1 .9|1 — (|1 . 10p . In this connection, we 
refer e.g. to [2], analyzing a dynamic model for frictional contact of a thermoviscoelastic body with 
a rigid obstacle, with the power law normal compliance condition for unilateral contact, and the 
corresponding generalization of Coulomb's law of dry friction. Unilateral contact is modelled by a 
normal compliance condition in [3] as well, where a dynamic contact problem for a thermoviscoelastic 
body, with frictional and wear effects on the contact surface, is investigated. A wide class of dynamic 
frictional contact problems in thermoelasticity and thermoviscoelasticity is also tackled in [24] . with 
contact rendered by means of a normal compliance law. 



«n < 0, ctn + Xun < 0, mn(ctn + Xmn) = in r c x (0, T), 



(1.14) 



<r T + Xu T | < c(d - ^ s )|ct n + Xm n |, 




(1.15) 
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In an extensive series of papers (cf. the references in the monograph [22]), C. Eck & J. 
Jarusek developed a different approach, which enabled them to prove existence results for dynamic 
contact problems, coupling dry friction and Signorini contact, without recurring to any regularizing 
operator. However, they used a different form of Signorini conditions, expressed not in terms of u 
but of dtu. They observed that this different law can be interpreted as a first-order approximation 
with respect to the time variable, realistic for a short time interval and for a vanishing initial gap 
between the body and the obstacle, and hence it is physically interesting. Within this modelling 
approach, existence results for contact problems with friction and thermal effects were for instance 
obtained in [33] and in [2"T] . 

Let us stress that, in all of these contributions the existence of solutions is proved for a 
weak formulation of the related PDE systems, which involves a variational inequality, and not an 
evolutionary differential inclusion, for the displacement. 

Analytical difficulties. From now on, we will take 6 = and thus confine ourselves to the 
investigation of the reversible case. We postpone the analysis of the irreversible case 5 > to 
the forthcoming paper |12j . where we will address also slightly different equations for $ and $ s 
obtained by a different choice of the bulk and surface entropy fluxes (here given by — V$ and — V$ s , 
respectively). For more details, we refer to the upcoming Remark 12.21 
Still, the analysis of system ([1.3111.12]) is fraught with obstacles: 

1) First of all, the highly nonlinear character of the equations, due to the presence of several 
singular and multivalued operators. In particular, since Coulomb friction is included in the 
model, a multivalued operator occurs even in the coupling terms between the equations for 
As, and u, and (|1.10|) features the product of two subdifferentials. 

2) A second analytical difficulty is given by the coupling of bulk and (contact) surface equations, 
which requires sufficient regularity of the bulk variables $ and u for their traces on T c to make 
sense. 

3) In turn, the mixed boundary conditions on u do not allow for elliptic regularity estimates which 
could significantly enhance the spatial regularity of u and <9tU. In particular, we are not in 
the position to obtain the H 2 (il; M. 3 ) -regularity for u and dtu. The same problem arises for 
$ , due to the third type boundary condition (|I.4[) . 

4) Additionally, the temperature equations (|1.3[) and (| 1 . 5|) need to be carefully handled because 
of the singular character of the terms dt ln(i?) and <9tln(?9 s ) therein occurring. In particular, 
they do not allow but for poor time-regularity of i? and i? s . 

Let us stress that the presence of several multivalued operators in system ([1.3111.12"]) is due to the 
fact that, all the constraints on the internal variables, as well as the unilateral contact conditions and 
the friction law, are rendered by means of subdifferential operators. Therefore, in our opinion the 
resulting formulation of the problem provides more complete information than those formulations 
based on variational inequalities. Indeed, it enables us to clearly identify the internal forces and 
reactions which derive by the enforcement of the physical constraints. 

In particular, let us dwell on the coupling of unilateral contact and the dry friction Coulomb 
law: as already mentioned, it introduces severe mathematical difficulties, unresolved even in the case 
without adhesion. These problems are mainly related to the fact that we cannot control ctn and <x t 
in ([1.9)1 and (| 1 . 10[) pointwise, due to the presence of two different non-smooth operators. For this 
reason, as done in [IT] and following Duvaut's work [20], we are led to regularize (| 1 . 10[) by resorting 
to a nonlocal version of the Coulomb law. The idea is to replace the nonlinearity in (|1.10p involving 
friction by the term 

c(^-i?,)|31(aJ(_ 00 ,o](tiN)|d(9 e u). (1.16) 
Accordingly, the coupling term c'($ — ■d s )dI(_ oc _ ](u^)\d t UT\ in (|1.4p and (|1.5[) is replaced by 

c'(tf - i9 s )|DI(9I ( _ 00 ,o ] (un))||^u t |. (1.17) 

In (|1.16p and (|1.17j) , Jl is a regularization operator, taking into account nonlocal interactions on the 
contact surface. We refer to [11] for further comments and references to several items in the literature 
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on frictional contact problems, where this regularization is adopted. In the forthcoming Remark 12.31 
we hint at the modeling derivation of the PDE system (jl. 3111.12"]) with the terms (|1.16|) and (|1.17p , 
while in Example 13 . 21 below we give the construction of an operator 51 complying with the conditions 
we will need to impose for our existence result. 

The difficulties attached to the coupling between thermal and frictional effects are apparent 
in the dependence of the friction coefficient c on the thermal gap $ — $ s . In order to deal with this, 
and to tackle the (passage to the limit in the approximation of the) terms (|1.16[) and (|1.17p . it is 
crucial to prove strong compactness for (the sequences approximating) $ and -d s in suitable spaces. 
A key step in this direction is the derivation of an estimate for d in BV(0, T; T4 Al,3+£ (fJ) / ) and for 
$ s in BV(0, T; W 1,2+t (T c )') for some e > 0. Combining this information with a suitable version of 
the Lions- Aubin compactness theorem generalized to BV spaces (see, e.g., [22]), we will deduce the 
desired compactness for (the sequences approximating) $ and i9 s . 

Our results. The main result of this paper (cf. Theorem [1] in Section l3~3|) states the existence of 
(global- in-time) solutions to the Cauchy problem for system 3111. 12]) . 

In order to prove it, we introduce a carefully tailored approximate problem for system (| 1 . 
I1.12[) . in which we replace some of the nonlinearities therein involved by suitable Moreau-Yosida-type 
rcgularizations. The existence of solutions for the approximate problem is obtained via a Schaudcr 
fixed point argument, yielding a local-in-time existence result. The latter is combined with a series of 
global a priori estimates, holding with constants independent of the regularization parameter, which 
enable us to extend the local solution to a global one. These very same estimates are then exploited 
in the passage to the limit in the approximate problem, together with techniques from maximal 
monotone operator theory and refined compactness results. In this way, we carry out the proof of 
Thm. [Q 

Plan of the paper. In Section [5] we outline the rigorous derivation of system 3111. 12")) and 

we comment on its thermodynamical consistency. Hence, in Section [31 after enlisting all of the 
assumptions on the nonlinearities featured in (|1.3II1.12"| and on the problem data, we set up the 
variational formulation of the problem and state our main existence result (cf. Theorem [1] below) . 
Section [4] is devoted to the proof of Thm. 14.11 ensuring the existence of (global-in-time) solutions to 
the approximate problem. Finally, in Section[5]we pass to the limit with the approximation parameter 
and conclude the proof of Thm. [TJ 

2 The model 

In this section we sketch the modeling approach underlying the derivation of our phase-transition 
type system for a contact problem, in which adhesion and friction are taken into account in a non- 
isothermal framework. The equations, written in the bulk domain and on the contact surface, are 
recovered by the general laws of Thermomechanics with energies and dissipation potentials defined 
in Q and on T c . They can be considered as balance equations, based on a generalization of the 
principle of virtual powers. Such a generalization accounts for micro-movements and micro-forces, 
which are responsible for the breaking of the adhesive bonds on the contact surface. We will not 
derive the PDE system (|1.3II1.12"[) in full detail, referring the reader to the discussions in [9] on the 
modeling of thermal effects in contact with adhesion, and in QT] for a (isothermal) model of adhesive 
contact with friction. Here, we will rather briefly outline the main ingredients of the derivation, and 
just focus on the modeling novelty of this paper. 

The phenomenon of adhesive contact is described by state and dissipativc variables, defining 
the equilibrium and the evolution of the system, respectively. The state variables are the symmetric 
strain e(u) , the trace of the vector of small displacements u on the contact surface, the surface 
adhesion parameter X , its gradient VX , and the absolute temperatures $ and i? s of the bulk domain 
and of the contact surface, respectively. The evolution is derived in terms of a pseudo-potential of 



5 



dissipation, depending on the dissipative variables df U, dtX, V#, Vi9 s , and the thermal gap on the 
contact surface I? — $ s . 



The free energy and the dissipation potential. The free energy of our system is written as 
follows 

* = * n + *r, 

"J n being defined in 51 and on T c . The bulk contribution "Jq is given by 

* S2 = 0(1 - ln(tf)) + 0tr(e(u)) + ij- (u)A' e e(u), (2.1) 

with Ag the elasticity tensor and tr(e(u)) the trace of e(u) . Note that, for the sake of simplicity, we 
have taken both the specific heat and the thermal expansion coefficient equal to 1 . The free energy 
on the contact surface is defined by 

*r =0a(l - + A(X)(0 a - ^ cq ) + I[o,i](X) + <t(X) (2-2) 

+ yX(w N ) 2 + yX|u T | 2 + /(_oo,0]("n) + y |VX| 2 , 

where cn, ct, k s are positive constants. Note that cn and Ct (which are the adhesive coefficients 
for the normal and tangential components, respectively) a priori may be different, due to possible 
anisotropy in the response of the material to stresses. However, for the sake of simplicity in what 
follows we let cn = ct = k s = 1 . In (|2.2[) , a is a sufficiently smooth function accounting for 
some internal properties of the adhesive substance on T c , such as cohesion: the simplest form for 
the cohesive contribution to the energy is u{X) — w s (l — X) for some w s > 0. Let us now briefly 
comment on the internal constraints. The energy is defined for any value of the state variables, but 
it is set equal to +oo if the variables assume values which are not physically consistent. Indeed, the 
indicator function iV—oo oi enforces the internal constraint un < , i.e. it renders the impenetrability 
condition between the body and the support. Finally, the term Im-a(X) forces X to take values in the 
interval [0, 1] . The function A provides the latent heat A' , while i? cq is a critical temperature, which 
governs the evolution of the cohesion of the adhesive substance with respect to the temperature. 

As far as dissipation is concerned, we consider in particular dissipative effects on the boundary 
due to friction. The pseudo-potential is 

$ = $ S2 + $ r 

$fi being defined in fJ and $r on r c . For the bulk contribution $n we have 

*n := \s{d t u)K v e{d t u) + ^|W| 2 , (2.3) 
while the contact surface contribution $r reads 

$r : = c(0 - 8 ) \-Rn + u N X\j(d t u) + hd t X\ 2 

(2.4) 

+ £/(-oo,o](d t X) + N^W| 2 + \k{X){d - tf s ) 2 , 

where the positive function c has the meaning of a friction coefficient, Rjq will be specified later (see 
(12.151) ). the function j is 

j(v) = |v T | for all v e K 3 . (2.5) 

5 > , and a and a s are the thermal diffusion coefficients in the bulk domain and on the contact 
surface, respectively. The positive (and sufficiently smooth) function k is also a contact surface 
thermal diffusion coefficient, accounting for the heat exchange between the body and the adhesive 
substance on T c . The assumptions on all of these functions have to guarantee that the pseudo- 
potential of dissipation is a no n- negative function, convex w.r.t. the dissipative variables (note that 
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•Q and i!} s are not dissipative variables) , and that dissipation is zero once the dissipative variables are 
equal to zero (see also the upcoming discussion on the thermodynamical consistency of the model) . 
Observe that, if 8 > 0, the model encompasses an irreversible evolution for the damage-type variable 
X, as it enforces dtX < 0. Furthermore, the 1 -homogeneity of the function j in (|2.5j) reflects the 
rate-independent character of frictional dissipation. 

Remark 2.1 Note that the functions a and A in the free energy (|2.2j) are related to the cohesion 
of the adhesive substance, as it results from Indeed, the term —a'(X) — X'(X)(-d s — $ oq ) is a 

(generalized) cohesion of the material, depending on the temperature. It represents a threshold for 
i|u| 2 to produce damage. 

The main novelty of this paper in comparison with is the fact that the friction coefficient 
c depends on the thermal gap i? — ?9 S (cf. [31] p. 12]). This relies on the modeling ansatz that friction 
generates heat (cf. |23p. so that ultimately we have the contribution c'(0 — $s)dI(-oo.o}( u ~N)\dt u T\ 
as a source of heat on the contact surface T c , both in the boundary condition (|1 .4[) for $ and in 
equation (|1.5[) for S . 



The balance equations. The equations for the evolution of the temperature variables are recov- 
ered from entropy balance equations in f2 and r c (cf. [15j for a detailed motivation of this approach 
and [9] for the application to contact problems), i.e. 

d t s + div Q = /). in £1 x (0, T), Q • n = F on 90 X (0, T) (2.6) 

with h an external volume source, F the entropy flux through the boundary, s the entropy, Q the 
entropy flux in the bulk domain, and 

d t s s + div Q s = F in T c x (0, T), Q s ■ n s = on dT c x (0, T), (2.7) 

where we have denoted by s s and Q s the entropy and the entropy flux on the contact surface, 
respectively. Note that, on r c the term F (involved in the boundary condition for (|2.6j) ) is the entropy 
provided by the adhesive substance to the body. We couple (|2.6p and (|2.7[) with the generalized 
momentum equations for macroscopic motions in Q and micro-movements in r c , viz. 

-div <x = f in n x (0,T), (2.8) 
ern = -R on T c x (0, T), u = on Y x x (0, T), ern = g on T 2 x (0, T) (2.9) 

where recall that er is the macroscopic stress tensor, f is a volume applied force, g is a known 
traction, — R is the action of the obstacle on the solid, and 

B - div H = in T c x (0, T), H • n s = on dT c x (0, T), (2.10) 

with H and B microscopic internal stresses, responsible for the damage of the adhesive bonds 
between the body and the support. 



The constitutive equations. To recover the PDE system p. 3111.12"]) . we have to combine (12.61) 
(|2.7[) with suitable constitutive relations, for the involved physical quantities, in terms of and <I>. 
Wc have 

9*n d* r 



<9i? ' 9?? s 

and 



(2.11) 



n--*?a o = _J^l. (212) 



As for as the flux through the boundary F , we impose 



9$r 

F = —, r on T r , F = on the remaning part. (2.13) 

9(0 -i9 a ) v ' 
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Then, we prescribe for (the dissipative and non-dissipative contributions to) the stress tensor 



nd, d OWQ OVfi 
<7 = (J + (J = 



de(u) de(d t u) ' 
and the reaction R = R^n + Rt reads 



Finally, B and H are given by 



B = B ,, + B .i = _^ + _j_ y (2 . 17) 



H H +H -av^ + 3(v(d t x))' (2 ' 18) 



(2.14) 



iijv = ^ = (2.15) 
Rt = R T d + R T = + ,^ $r . . (2.16) 



Remark 2.2 Let us briefly comment on the entropy flux laws for Q and Q s : as shown by (|2.3[) 
and (|2.4p , they depend on the choice of the thermal diffusion coefficients a and a s . In general the 
latter are functions of the temperature variables -d and i? s (cf. e.g. [22]). If a(d) = a > we get 
Q = -aViJ. In this case, for the heat flux q = #Q we obtain q = — a$V?9 = — ^Vt? 2 . Analogously, 
if a($) = 4 we get Q = — Vln(i?) and for q the standard Fourier law q = — Vi?. Hence, in the case, 
e.g., a($) = $ (which is admissible for thermodynamics) we have Q = — $Vi? = — ^Vi9 2 . Analogous 
considerations hold for Q s . 

In what follows, we fix as heat flux laws 

a(0) = a=l, a s (i? g ) =a g = 1. (2.19) 

Furthermore, as already mentioned, we will confine our analysis to the case of a reversible evolution 
for the damage-type variable X , taking 

5 = in 3> r . 

With these choices, combining (|2TTT]) - ([2TT8]| with (j2~6| - ([2~T0)) we derive the PDE system (p^3l [TTT2]) . 

Remark 2.3 Let us point out that, the PDE system tackled in this paper, featuring the non-local 
rcgularization for the Coulomb law 51, can be derived following the very same procedure described 
above. Indeed, it is sufficient to replace $r by 

$r : = c(0- 0.) \K(-Rn + u N X)\j{d t u) + ^\d t X\ 2 

a (4) 1 (2 ' 20) 

+ «(-oo,o](3 t X) + ^^|Vtf s | 2 + - S ) 2 . 

Thermodynamical consistency. Let us now briefly comment on the derivation of the model in 
relation to its thermodynamical consistency. In particular we discuss the validity of the Clausius- 
Duhem inequality for the whole thermomechanical system we are dealing with. The latter consists of 
the body f2 and the contact adhesive surface r c , encompassing the dissipative heat exchange between 
the body and the adhesive substance on T c , due to the difference of the temperatures. Recall that 
we have derived the evolution equations for the temperature variables from an entropy balance in 
the domain D, and on r c . To prove thermodynamical consistency (which corresponds to showing 
that dissipation is positive), we have to discuss the entropy balance for the whole system, i.e. also 
including the interactions between the bulk domain and the adhesive substance. 
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As far as the domain £7 is concerned, this property is fairly standard (the reader may refer to 
|26j): it is ensured once $n has the characteristics of a pseudo-potential of dissipation. In particular, 
this requires that, in (|2.3p the tensor K v is positive definite, and the thermal diffusion coefficient a 
is non-negative. This leads to positivity of the dissipation, viz. 

e(d t u)K v e(d t u) +a(i?)|Vi9| 2 > 0. 

Then, the Clausius-Duhem inequality is obtained by writing the energy balance and exploiting the 
already specified constitutive relations. 

Second, the internal energy balance on T c is recovered by the first principle, in which the 
effective power of the internal forces responsible for the damage process is included. Again using the 
constitutive relations, it yields 

s d t s s + s div Q s = 6 S Q S + B d d t X + U d ■ \7d t X - Q s ■ W s , (2.21) 

where we denote by Q s the entropy received by the adhesive substance from the solid and B d , H d , Q s 
are defined by (|2.17|) . (|2.18[) . and (|2.12[) . Note that Q s actually represents an entropy source. In our 
model we have prescribed that (see (|2.6[) - (|2.7[1 ) 

-Q = Q S = F, (2.22) 

Q denoting the entropy received by the solid from the adhesive substance. Note that by (|2.6[) we 
have 

Q-n + Q = 0. (2.23) 

To deduce the conditions on F which ensure the thcrmodynamical consistency, we write an 
entropy equality on the contact surface, accounting for the interaction between the body and the 
adhesive substance, i.e. 

ms in t = -Q#-Q s -&s, (2.24) 

where Si nt is the entropy exchange, defined in terms of the variable d — i(?9 + $ s ) (for any further 
details sec [26]). Now, the Clausius-Duhem inequality on the contact surface (combining the adhesive 
substance and its interaction with the body) reads 

d t s s + div Q s + dtSint > Q • n. (2.25) 

Thus, exploiting (pOTj) . (|2~22"j) . (f2T23|) . and ([2~24)) . inequality (|2~25|) follows once it is ensured that 
(recall that the absolute temperatures are strictly positive) 

F{6 - 6 S ) + B d d t X + H d • Vd t X - Q s • V6» s > 0. 

Observe that this corresponds to prescribing in (|2.4[) that 

a s >0, c>0, fc>0, c'(i? -i? s )(i9- S ) >0. 

Indeed, all of the above conditions are ensured by (|2.19j) and by the forthcoming Hypotheses (I) and 
(V) on c and k. 

Remark 2.4 Note that, in the case no heat is exchanged between the body and the adhesive sub- 
stance, the Clausius-Duhem inequality on the contact surface is granted upon requiring that 

B d d t X + H d • Vd t X - Q s ■ V6» s > 0. 

This follows from (|2"2T|) taking Q s = . 
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3 Main results 



3.1 Setup 

Throughout the paper we shall assume that 

57 is a bounded Lipschitz domain in R 3 , with 
dil — Ti U T2 U r c , Ti, i = 1, 2, c, open disjoint subsets in the relative topology of dil, such that 
M 2 (ri), J{ 2 (r c ) > 0, and r c cl 2 a sufficiently smooth flat domain. 

(3.1) 

More precisely, by flat we mean that T c is a subset of a hypcrplane of R 2 and !K 2 (r c ) = L 2 (T C ), 
L d and "K d denoting the d -dimensional Lcbesguc and Hausdorff measures. As for smoothness, we 
require that r c has a C 2 -boundary: thanks to, e.g., [27l Thm. 8.12], this will justify the elliptic 
regularity estimates we will perform on the solution component X . 

Notation 3.1 Given a Banach space X , we denote by (■, -) x the duality pairing between its dual 
space X' and X itself and by || ■ ||x the norm in X. In particular, we shall use the following 
short-hand notation for function spaces 

H:=L 2 (Cl), V 1=^(0.), H := L 2 (f2;R 3 ), V := if x (f2; R 3 ), 
H Tc := L 2 (r c ), Vr c := H 1 ^), Fr c := ffjj^ ( r c)> 
W := {v 6 V : v = a.c. on TJ, H Fc := £ 2 (r c ;R 3 ), Y Fc := ^ ri (r c ;R 8 ), 
where we recall that 

F oo,r, ( r <0 = { w e -^ 1/2 (r c ) : 3ti£ H 1/2 {T) with w = w in T c , w = in Ti | 

and Hqq 2 Ti (r c ; R 3 ) is analogously defined. We will also use the space Hqq 2 Ti (T2; R 3 ). The space W 
is endowed with the natural norm induced by V . We will make use of the operator 

A:Vr c ->Vr c (AX,w) Vr := J VXViudx for all X, w G V Ta (3.2) 
and of the notation 

m(w) := £^r^\ y wdi for u> G L 1 (A). (3-3) 

Useful inequalities. We are going to exploit the following trace results and continuous embeddings 
VcL 4 (r c ), WcL 4 (r c ;R 3 ), Vr € C L q (T c ) for all 1 < q < 00, (3.4) 
and the fact that, by Poincare's inequality, for every Lipschitz domain A C R d 

3C>0V«6^); \\v\\w^(A)<C(\\v\\ LHA) + \\Vv\\ LHA) ), (3.5) 
where H|z,i(A) can be replaced by |m(u)|. 

Linear viscoelasticity. We recall the definition of the standard bilinear forms of linear viscoelastic- 
ity, which are involved in the variational formulation of equation (|1.7p . Dealing with an anisotropic 
and inhomogeneous material, we assume that the fourth-order tensors K e = (a.ijkh) and K v = 
(bijkh) , denoting the elasticity and the viscosity tensor, respectively, satisfy the classical symmetry 
and ellipticity conditions 

o>ijkh = a-jikh = a-khij , bijkh = bjikh = bkhij for k, h = 1, 2, 3 
3 a > : aijkhiij^kh > "oCy6j V&j : = for i,j = 1, 2,3 , 
3 (3 > : bijkh£ij£kh > Po&j&j V&j : = for i, j = 1, 2, 3 , 
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where the usual summation convention is used. Moreover, we require 

Oijkh,kjkh £ L°°(Q,), i,j,k,h= 1,2,3. 

By the previous assumptions on the elasticity and viscosity coefficients, the following bilinear forms 
a, b : W x W -> E, defined by 

a(u,v) := / Oijfcfcefcfc(u)e^(v)da; forallu,veW, 

6(u,v) := / bijkhSkhiu^ijirfdx for all u, v 6 W, 

turn out to be continuous and symmetric. In particular, we have 

3C>0: |a(u,v)| + |6(u,v)| < <7||u|| w ||v||w forallu,veW. (3.6) 

Moreover, since Ti has positive measure, by Korn's inequality we deduce that a(-, •) and &(•, •) are 
W -elliptic, i.e., there exist C a , Cb > such that 

o(u,u) > CalluH^v forallueW, (3.7) 
6(u, u) > C b ||u||w for all u e W. (3.8) 

3.2 Assumptions 

We specify all of the assumptions on the nonlinearities in system (jl. 3111.121) . 
Hypothesis (I). As for the friction coefficient c : R — > (0, +oo) , we require that 

ceC^M), 3ci,c 2 >0ViGl : c(i)>d, |c'(x)| < c 2 , c'(x)x>0. (3.9) 

For instance, the function c(x) = J Q arctan(i) dt + c\ = a;arctan(a;) — 2"^ n (^ x ^ Cl corn P ncs 
with (|33|) . 

Hypothesis (II). We generalize the operator -00,01 by replacing It-cofi] m (|2-2p with a function 

</> : R — > [0, +00] proper, convex and lower semicontinuous, with <fi(0) = (3.10) 

and effective domain dom(0) ; let us emphasize that the physical case, in which the constraint mn < 
on T c is enforced, occurs when dom(</>) C (— 00, 0] (and it is included in our analysis). Then, we 
define 

<P ■ Yr c -> [0, +00] by ip{v) := i J1 ° (3.11) 

I +00 otherwise. 

Hence, we introduce 

if : Yr c [0, +00], defined by <p(u) := v?(u N ) for all u e Y Fc . (3-12) 

Since y : Yp c — >• [0, +00] is a proper, convex and lower semicontinuous functional on Yp c , its 
subdiffcrcntial dip : Yr c =4 Y r is a maximal monotone operator. 

Hypothesis (III). Concerning the regularizing operator ft, in the lines of [TT] we require that there 
exists v > such that 

ft : L 2 (0,T; Y r J L°°(0, T; L 2+ "(r c ; R 3 )) is weakly-strongly continuous, viz. 
77 n -T7 inL 2 (0,T;Y rc ) ft(rj n ) -> in L°°(0, T; L 2+,y (r c ; R 3 ) 

for all (rj„), rj € £ 2 (0, T; Y r J . It is not difficult to check that (pTT3|) implies that ft : L 2 (0, T; Y r J -> 
L°°(0, T; L 2+y {T c ; R 3 )) is bounded. In Example 1 3. 2 1 at the end of this section we are going to exhibit 
an operator ft complying with (|3.13l) . 
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Hypothesis (IV). We generalize 9/[o,i] in (jl.lip by considering the multivalued operator d/3 : 
[0, +00) =4 R, with 

B : R — > (— 00, +00I proper, convex and lower semicontinuous, 

(H4) 

such that dom(/3) C [0,+oo). 
In what follows, we use the notation B := 38. 

Hypothesis (V). We assume that the nonlinearities k in Ql.4p - p.5p . A in Q1.5P and (jl.lip . and a 
in (|LTT|) fulfill 

k : R — > [0, +00) is Lipschitz continuous, (3-14) 
A 6 (^(R), with A' : R ->• R Lipschitz continuous, (3.15) 
a G C^R), with cr' : R -> R Lipschitz continuous. (3.16) 

Assumptions on the problem and on the initial data. We suppose that 

h G L 2 (0,T;V') nL\0,T;H) , (3.17) 
feL 2 (0,T;W), (3.18) 
geL^OJ;^^;! 3 )'). (3.19) 

It follows from ([3TT5|1 - (|339]) that, the function F : (0, T) -> W defined by 

(F(i),v) w := (f(f),v) w + (g(i),v) 1/2 rr , H3 s for all v G W and almost all t G (0,T), 

satisfies 

FeI 2 (0,T;W). (3.20) 
Finally, we require that the initial data fulfill 

tfoeL 1 ^) and Ln(i?o) G # , (3.21) 

^eL 1 ^) and ln(^°)Gi?r c , (3.22) 

uo G W and uo G dom(<p) , (3.23) 

X eVr c , ^(X )Gi 1 (L c ). (3.24) 

We conclude with the example, partially mutuated from |111 Ex. 2.4], of an operator 31 complying 
with Hypothesis (III). 

Example 3.2 Fix £ : T c x T c -> R 3 such that i G L 2+V (T C ; Y r J for some z/ > 0, and for all 
T7 G L 2 (0,T; Y' r J set 

3?(T7)(x,t) := ^jf ( V {;s),e(x,-))v re ds^w for a.a. (x,t) G T c x (0,T), 

where w is a fixed vector, e.g. with unit norm, in R 3 . Then, for almost all (x, t) G T c x (0, T) there 
holds 



|IR(T7)(x,i)| < / |(r7(., S ),^(x,-)} Yr Jd S <t 1 / 2 |K(x,-)|| Yr J|r 7 |U2 (0 , t:Y , (3.25) 
Jo 

Integrating (|3.25p over T c , we easily conclude that for all t G (0,T) 

\\%(v)(;t)\\L* + »(T c ;R3) <T 1/2 ||»?l|w(o I r ; Yt e )ll^ll£»+-'(r ei Y r .) ) 

hence 31 : L 2 (0,T;Y^) -> L°°(0, T; L 2+l/ (r c ; R 3 )) is a linear and bounded operator. Furthermore, 
it fulfills lpTT3]) . Indeed, let rj n T7 in L 2 (0,T;Y^ ): for almost all (x,t) G T c x (0,T) we have 



R(Vn)( x , t )= ( J 7n(-> S ) I l(o,*)^(a;. -)}Y rc w ds -> / faO, s), l(o,t)^(a;, -))Y rc w ds = £(77) (x, i) 
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as n — > oo. Then, estimate (|3.25p and the dominated convergence theorem yield 3l(r) n ) —> !R(tj) 
in L9(0,T;L 2+,y (r c ;R 3 )) for all 1 < q < oo. In order to conclude that 3t(r] n ) -> ft(fj) in 
L°°(0, T; L 2+ly (r c ; M 3 )) , it is sufficient to observe that the sequence {R{T) n )) n is in fact compact 
in C°([0, T]; L 2+lJ {T c ] R 3 )) . This follows from the Ascoli-Arzela theorem, since, arguing as for (|3.25[> . 
it is not difficult to see that, if (f] n ) n C L 2 (0, T; ~Y' r ) is bounded, then (jR(rj n )) n fulfills the equicon- 
tinuity condition for all < s < t < T 

p(r,J(.,t) - X(Vn)(; s)IU^(r c;R 3) < (t - s) 1/2 ||77 JU=(,, Wo )IKIU^(r o;Yrc ) < C(t - s) 1 ' 2 . 

3.3 Statement of the main result 

We now specify the variational formulation of system p. 31 11.121) . 

Problem 3.3 Given a quadruple of initial data t9°, Uo, Xo) fulfilling (|3.2ip - (|3.24j) . find a seven- 
tuple ($, $ s ,u, X, 77, with 

G i 2 (0, T; V) n L°°(0, T; L 1 (f2)) , (3.26) 
Ln(i?) G L co {0,T;H)r\H 1 (0,T;V') , (3.27) 
S G L 2 (0, T; F r J n L°°{0, T; L\T C )) , (3.28) 

in(0») Gi 00 (o I r;fr rB )nfr 1 (o J r;V r / e ), (3.29) 

ueff'^TjW), (3.30) 

X G L 2 (0, T; H 2 {T C )) PI L°°(0, T; Vr B ) n T; #r c ) , (3.31) 

f|GL 2 (0,r;Yt e ), (3.32) 

/x G L 2 (0,T;HrJ, (3.33) 

£eL 2 (0,T;ff r J, (3.34) 

satisfying the initial conditions 

0(0) = i9 a.e. in fi , (3.35) 

tf s (0) = tf° a.e. inT c , (3.36) 

u(0) = u a.e. in il , (3.37) 

X(0)=X o a.e. in T c , (3.38) 



and 



(d t ln('&),v) v - / div(<9 t u)t>d:r + / V0Vvdx + / fc(X)(tf - d s )v dx 
Jn Jn Jr c 

+ [ c'(i?-i? a )|3l(i7)||5 t u T |t;da;= (h,v) v VweV a.e. in (0,T) 
Jr c 

(d t \n(d s ),v) Vr - / d t X(X)vdx+ [ \7d s Vvdx 
Ta Jr c Jr c 



(3.39) 



(3.40) 



(3.41) 



k(X)($ -$ s )vdx + [ c'(i?-i? g )|3l(77)||afUT|uda; \/veVr c a.e. in(0,T), 
Jr c 

b(dtu, v) + a(u, v) + / 0div(v)d;r + / Xu-vdx 
Jn Jr a 

+ (T7,v) Yr + / c(-&-'& a )fi-vdx = (F,v) w forallveW a.e. in (0,T), 
Jr. 

77 e cV(u) in Y^, a.e. in (0,T), (3.42) 

= |ft(j?)|z with z G d(d t u) a.e. in T c x (0,T), (3.43) 

d t X + AX + t + a'{X) = -\'{X)d s - i|u| 2 a.e. in T c x (0, T) , (3.44) 

£ G /3(X) a.e. in T c x (0,T). (3.45) 
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Note that, to simplify notation we have incorporated the contribution — A'(X)$ oq occurring in 
into the term a'(X) in ([3T4|) . 

Remark 3.4 Observe that the subdifferentials in (|3.43p and p.45|) are multivalued operators on 
R 3 with values in R 3 , and on R with values in R , respectively. Hence the related subdifferential 
inclusions for z and £ hold a.e. in T c x (0, T) . Instead, (|3.42[) features the (abstract) operator 
dip : Yr r =4 Y Fc , thus (|3.42[) holds in Yp c . Note that, in the case wc further assume (for physical 
consistency) that dom0 C (— oo,0], it still follows from the definition Q3.12p of cp that, if u G 
dom(<9y), then u complies with the constraint un < a.e. in T c x (0,T). 

We are now in the position to state our existence theorem for Problem 13.31 Observe that we obtain 
enhanced regularity in time for i) and ?? s thanks to some refined BV -estimates (cf. Remark 14.101 
later on). Relying on (|3.13[) in Hypothesis (III), we also find (|3.48[) for /j, . 

Theorem 1 

In the framework of (|3.ip . under Hypotheses (I)-(V) and conditions (|3 . 1 7[) - (|3.19p on the data h, f, 
g, and (|3.2ip - (|3.24p on t9 , §° s , u , and Xo , Problem [373[ admits at least a solution (<&, $ s , u, X, r/, fJ,,£), 
which in addition satisfies 

d G BV(0, T; W hq {ny) for every q>3, (3.46) 
i? a G BV(0, T; W 1,ff (r c )') for every a > 2, (3.47) 
/j,eL° o (0,T;L 2 + u (T c ;W s )) with v > from (TO]) . (3.48) 

The proof will be developed throughout Sections 0H5] First, we will analyze a suitable approximation 
of Problem 13.31 for which we will obtain the existence of solutions in Thm. [431 In Section [5] we will 
then pass to the limit and show that (up to a subsequence) the approximate solutions converge to a 
solution of Problem 13. 31 

Notation 3.5 In what follows, we will denote most of the positive constants occurring in the calcu- 
lations by the symbols c, c', C, C , whose meaning may vary even within the same line. Furthermore, 
the symbols 7j, i = 0, 1, . . . , will be used as place-holders for several integral terms popping in the 
various estimates: we warn the reader that we will not be self-consistent with the numbering, so that, 
for instance, the symbol I\ will appear several times with different meanings. 

4 Approximation 

Here we focus on the approximation of the PDE system (|3.39l I3.45P through Yosida-type regulariza- 
tion of some of the (maximal monotone) nonlinearities therein involved. For the related definitions 
and results we refer to the classical monographs |SJ [TB] ■ 

In Sec. 14.11 we justify the regularizations we will perform, giving raise to the approximate 
Problem (P e ) . The existence of solutions is proved in two steps: first, in Sec. 14.21 we show that 
Problem (P e ) admits local- in-time solutions by means of a Schauder fixed point argument. We then 
extend them to solutions on the whole [0, T] relying on the global a priori estimates which we derive 
in Sec. 14.31 In fact, we will obtain estimates independent of the approximation parameter e: they 
will be the starting point for the passage to the limit as e J, developed in Section \5\ 

4.1 The approximate problem 

In order to motivate the way we are going to approximate Problem 13 - 31 let us discuss in advance 
some of the global a priori estimates to be performed on system (|3.39l - l3~45|) . Clearly, these estimates 
correspond to (some of) the summability and regularity properties (|3.26p - (|3.34p required for solutions 
to Problem [3?3l As we will see, the related calculations can be developed on system (|3.39ll3T5l) only 
on a formal level: they can be made rigorous by means of the Yosida-type regularizations we will 
consider (cf. the global a priori estimates on the approximate solutions performed in Section 14.31) . 
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Heuristics for the approximate problem. The basic energy estimate for system (|3.39ll3~TT)l) 

consists in testing (jcT3l)j) by ■&, p~40j) by -d s , ([3~4T|) by d t u, (|X44| by d t X, adding the resulting 
relations, and integrating in time. Taking into account the formal identities 

/ {dtH#),*)v dr = ll*(*)lk(n) ~ WMlhu) , 

J \ (4-1) 

(dt H# s ),ti s ) Vrc dr = |K(i)|Ui ( r c ) - iKlkMr.) , 

this estimate leads, among others, to a bound for z? in L°°(Q, T; L 1 (fi)) and, correspondingly, for -d s 
in L oo (0,T;L 1 (r c )) (cf. (|3~26| and (|3~28| ). As a first step towards making (j4~Tj) rigorous, following 
[H [10] we will 

1. replace the logarithm In in equations p.39[) and (|3.40p by its approximation L £ : R — > M 

L £ {r) :=er + \n £ (r), (4.2) 

where for e > ln e denotes the Yosida regularization of the logarithm In. Therefore, L £ is 
differentiable, strictly increasing and Lipschitz continuous, see also Lemma 14.21 below. 

As pointed out in Remark 14.121 below, this procedure is not sufficient to justify (|4.1[) completely. 
In order to do so, following JSJ [TO] we should also add a viscosity term both in (|3.39p and (|3.40p , 
modulated by a second parameter v . We choose to overlook this point for the sake of simplicity. 
Anyhow, let us stress that approximating the logarithm In by L £ makes the test of (|3.39[) by ln($) 
rigorous (and, respectively, it justifies the test of (|3.40l) by ln(^ s )). This gives raise to an estimate 
for ln(tf) in L°°(0,T;H) (for ln(?? s ) in L°°(0, T; H r J , respectively). 

A consequence of the aforementioned energy estimate and of a comparison argument in the 
momentum equation p.4ip is the following estimate 

\\c(#-4,)n + ri\\ L n 0i T^ e )<C. (4.3) 

Using the formal observation that \x and r) are orthogonal, from (|4.3p one concludes the crucial 
information 

||c(??-^ s )/x|| L 2 (0>x . no) + ||77|| L 2 (0>x . Y /,j <C. (4.4) 

In order to justify this argument, we need to suitably approximate the maximal monotone operator 
dip : Yr c Yp c in such a way as to replace r] s Yp c in p.4ip with a term T] e having null 
tangential component, cf. (|4.5|) below. Following |llj . we will thus 

2. replace the function </>, which enters in the definition of the functional ip through (|3.11l) 
and (|3.12p . by its Yosida approximation <fr £ : R —> [0, +oo) . We recall that (f> e is convex, 
differentiable, and such that <f>' E is the Yosida regularization of the subdiffcrential d<f> : R R . 

Therefore, in this way we will replace the constraint (|3.42p by its regularized version 

r] e := <t>' e (u N )n a.e. in T c x (0,T). (4.5) 
Furthermore, we introduce the function J £ : R — > R defined by 

3 £ (x) := f sL' £ (s)ds (4.6) 
Jo 

which we will use in the proof of the existence of solutions to Problem (P E ) , and in particular in 
the derivation of suitable a priori estimates. In the two following lemmas, we collect some useful 
properties of the functions ln e , L £ , and J £ . 
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Lemma 4.1 The following inequalities hold: 



3e»>0: Ve e (0,e») Vi > In'(x)<-. 

x 

Ve > VxeR hi' e (x) > 



As a consequence, the junction 3 S 



\x\+2- 
(|4.6p satisfies 



3e*>0: Vee(0,e*) Va; > !L(V) < -a; 2 + 2x ; 
3Ci,C 2 >0: Ve>0 Viet %{x) > |x 2 + d|ar| - C 2 



(4.7a) 
(4.7b) 

(4.8a) 
(4.8b) 



Proof. Estimates (|4.7a|) and (|4.7b[) have been derived in [TUl Lemma 4.1] (cf. also [331 Lemma 4.2]). 
Inequalities (|4.8a[) and (|4.8L>|) can be deduced from (|4.7a|) — (|4.7b[) by arguing in a completely analogous 
way as in the proof of [101 Lemma 4.1], to which the reader is referred. ■ 



Lemma 4.2 The function L £ 



satisfies: 



e < L' e {x) < s + - for all x e 
1 1 



L' E (x) L' s {y) 



^ \ x ~ y\ f or a M x i y g 



(4.9) 
(4.10) 



Proof. The left-hand side inequality in (|4.9p directly follows from the definition (|4.2[) of L £ and 
from estimate (|4.7b>[) for ln^. . In order to prove the right-hand side estimate, it is sufficient to show 
that 

(4.11) 



(4.12) 



\L e (x) -£>e(y)\ < I e + -J \x - y\ for all x,y e K. 
This can be checked by recalling that the Yosida-rcgularization ln e is defined by 

\n e (x) := -(x - p e (x)) 



where p £ = (Id+eln) -1 : K — > (0,+oo) is the £ -resolvent of In. Plugging (|4.12[) into the definition 
(|4.2[) of L e and using the well- known fact that p e is a contraction, wc immediately deduce (|4.11l) . 

Observe that, by the first of (I4.9p the function x t-> r ,\ ; is well-defined on M. In order to 
show that it is itself a contraction, we use the formula 



\n e (x) 



1 



p e {x) + £ 

(cf. p~3l Lemma 4.2]). Therefore, for every x, y G ffi 



1 



1 



p e (x) +£ 



for all x £ 



Pe(y) + £ 



(4.13) 



z 2 + ep e {x) + l £ 2 +ep £ (y) + l 

\Pe(x)-p e (y)\ 



< \p e (x)-p e (y)\ 



(e 2 + ep e (x) + l)(e 2 + e Pe {y) + 1) 
and (|4.10j) ensues, taking into account that p e is a contraction. ■ 

Furthermore, we will supplement Problem (P £ ) by the approximate initial data (i?q, i9® ,£ , Uo, Xq) , 
where the family (^g,?9° ,e ) e fulfills 



(0§) e C ff, «' e ) e C iJ r 



(4.14a) 
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and it approximates the data ($0,$°) from (|3.2ip - (|3.22l) in the sense that 

(0§, i?°' e ) -> (i?o, #°) in x Ll ( r c) as £ I , (4.14b) 

||ln e (^)|| H <||m(tf )lk, ||ln e (^ £ )|| ffrc < ||ln(^)|| Hrc foraU £ >0, (4.14c) 
3S>0: Ve>0 e(||^|k + H^lkrJ < (4.14d) 

3S!>0: Ve>0 /" J e (0§(s)) dx < + ||0o|Ui ( n)) . (4.14e) 
3S 2 >0:Ve>0 jf a £ (^ £ )dx<S 2 (l+|[i?°|| il(rc) ). (4.14f) 

Indeed, (|4.14a|) reflects the enhanced regularity (|4. 1 9[) and (|4.21[) required of solutions (#, # g ) to 
Problem (P E ) . In what follows, we give an example of construction of a sequence (i^g, i?° ,6 ) e fulfilling 
properties (|4.14[) . 

Example 4.3 We will carry out the construction of the sequence (i?§)e only, the argument for (■d® ,s ) e 
being completely analogous. For all e > and a.e. in , let us define 

i9q := min-j^o,^"} for some a > . (4-15) 

Observe that > a.e. in (being i? > a.e. in Q thanks to the second of (|3.2ip ). tf§ e L°°(Q) , 
and moreover 

Hence, (|4.14b[) is an immediate consequence of the dominated convergence theorem. Next, noting 
that H^olllf — 1^1 £~ 2q and choosing a = i , we can deduce 

3S>0: Ve>0 e 1/2 \\tf E \\ H < S (4.17) 

and also (|4.14d|) follows. Moreover, ()4.8a|) . (|4.16|) . and (|4.17|) give (|4.14e|l . Finally, relying on the 
definition (|4.15[) and on well-known properties of the Yosida regularization ln £ , we find 

| ln e (^)| < | ln £ (^ )| < I M^o)| a.e. in Q (4.18) 

whence (|4.14c|) . 

All in all, the approximation of Problem 13. 31 reads: 

Problem 4.4 (P E ) Let us a consider a quadruple of initial data (i9g, ??°' e , Uq, Xq) satisfying (|3.23p - 



(j3~24l) and (|4~14j) . Find a sextuple (tf, tf s , u, X, /*,£), fulfilling 

•d e L 2 (0,T;V) nC°([0,T];F), (4.19) 

JC e (#) Gi 2 (0,T;F)nC ([0,r];iJ)ni/ 1 (0,T;F') I (4-20) 

d s e L 2 (0,T;V r J n C°([0,T}; HrJ, (4.21) 

L £ (# s ) e i 2 (0, T; VrJ D C°([0, T]; F r J D ff X (0, T; t£J, (4.22) 



and such that (u, X, £t,£) comply with (|3.30|) - (|3.31|l . Q3.48p . (|3.34[) . satisfy the initial conditions 

i?(O) = 0g a.e. in f2, (4.23) 
tf s (0) = ??°> e a.e. inr c , (4.24) 
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as well as (|3.37p - p.38p . and the equations 

{dtCe{$),v) v - I dw{d t \i)vdx+ I V&Vvdx+ f k(X)(ti - fi s )v dx 

+ / c / (i?-t? a )|aJ(^(ti N )n)||8 t u T |ud2;= (h,v) v V-u G V a.e. in (0,T), 



(d t L £ (d a ),v) Vrc - d t X{X)vdx + V$ s Vvdx 
Jr c Jr. 

k(X)(tf -tf s )vdx + c'(i9 - tf s )|ft(<^(w N )n)||<9 t u T |vda; Vv 6 Vr c a.e. in(0,T), 



6(9tU, v) + a(u, v) + / ^div(v)dx+ / Xu-vdx 

° c (4.27) 

+ y ^(u N )n'vds+ y c(i?-i? s )/i'vdi= (F,v) w for all v G W a.e. in (0,T), 

/x = |ft(^(u N )n)|z with z 6 d(<9 t u) a.e. in T c x (0,T), (4.28) 

as well as relations (|3.44[) — (|3.45[) . 

Observe that, in (|4.25|) , (|4.26|) , and (|4.28p the term ft(<^(u N )n) needs to be understood as ft( J(^(u N )n 
where J denotes the embedding operator from L 2 (0, T; L 4 (r c ; M 3 )) to L 2 (0, T; Y r J . 

In the next two sections, we address the existence of solutions to Problem (P e ) for e > 
fixed. That is why, for notational convenience we will not specify the dependence of such solutions 
on the parameter e . 

4.2 The approximate problem: local existence 

The main result of this section is the forthcoming Proposition 14.91 stating the existence of a local- 
in-time solution to Problem (P e ) . We prove it by means of a Schauder fixed point argument, which 
relies on intermediate well-posedness results for the single equations in Problem (P £ ) . 

Fixed point setup. In view of Hypothesis (III), we may choose S € (0, 1) such that 

ft : L 2 (0, T; Y r J -t L°°{0, T; L^J (r c ; M 3 )) is wcakly-strongly continuous (4.29) 
(and therefore bounded). For a fixed r > and a fixed constant M > 0, we consider the set 
y r = {(0, ti s , u, X) G L 2 (0, r; iT 1 " 5 ^)) x L 2 (0, r; H 1 " 5 ^)) x L 2 (0, r; ff 1 - 5 ^; R 3 )) x X T : 

II^II^CO^;^ 1 - 5 ^)) + ||'#s]].L2(0,-nif 1 -' 5 (r c )) + ||u||i2(0, T ;i3-l-«(n ; ]R3)) + ||X|| L 2 (0:T;fl - r j < M, 

with the topology of L 2 (0,t; H 1 -* (ty) x L 2 (0, r; H 1 ' 8 ^^) x L 2 (0, r; H 1_<5 (fi; 1R 3 )) x L 2 (0, r; ff r J , 
where 

I r :={XeL 2 (0,T;ff r J : X G dom(/3) a.e. on T c X (0, r)}. 

We are going to construct an operator T mapping ^ into itself for a suitable time < T < T , in 
such a way that any fixed point of T yields a solution to Problem (P e ) on the interval (0, T) . In the 
proof of Proposition 14. 9[ we will then proceed to show that T : — s- does admit a fixed point. 

Notation 4.5 In the following lines, we will denote by Si, i = 1,2,3, a positive constant depending 
on the problem data, on Si and S 2 in (|4. 14e)) - (|4. 14f|) . on M > in (|4.30p . but independent of e > 0, 
and by Si(e) a constant depending on the above quantities and on e > as well. Furthermore, 
with the symbols Wi(A) , Ki j(A), ... , we will denote the projection of a set A on its i-, or 
componcnt, .... 
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Step 1 : As a first step in the construction of T, we fix -d s , u, X) € y T and prove a well-posedness 
result for (the Cauchy problem for) the PDE system (|4^7ti4~28l) . with (&, d s , X) in place of (t?, t?„, X) , 
and 

c(?-5.)|0l(^(«N)n)| replacing c(0 - 0,)|:R(#(u N )n)|. 



Lemma 4.6 Assume (|3.ip . Hypotheses (I)-(III), and suppose that f, g, Uo comply with (|3.18|) - 
([3TT9]) . and (|3~23]) . 

Then, there exists a constant S± > sitc/i that for all (??, i? s , u, X) G ^ T iftere exists a unique 
pair (u, /x) G ff^TjW) x L°°(0, r; L 2+,y (r c ; M 3 )) fulfilling the initial condition (|3~3"?|) and 



b(dtu, v) + a(u, v) 



i? div(v) dx 



Xu ■ v dx 



<^(u N )n-vdx + / c(tf-?? s )/Li- vdz = (F,v) w for all w a.e. in (0,T), 

= |3i(0^(uN)n)|z z £ d(<9 f u) a.e. in T c x (0, T), 



(4.31) 



and £/ie estimate 



lullffHO^W) + ll z llL°°(r c x(0,r)) + II A* llL oo (0,T;L 2 + "(r c ;R 3 )) < -Si. 



(4.32) 



Proof. Observe that (|4.31[) has the very same structure of the momentum equation in the isothermal 
model for adhesive contact with friction tackled in [llj , for which the existence of solutions was proved 
by passing to the limit in a time-discretization scheme. The arguments from [llj can be easily adapted 
to the present setting, also taking into account that the term 



01 := c(d -d s )\Jl((f>' E {u N )n)\ is in L 2 (0, r; L 4/3 (r c )). 



(4.33) 



Indeed, P~33"]) can be checked by observing that c(&-d s ) G L 2 (0, t; L A / {l+2S \Y c )) by trace theorems 
and the Lipschitz continuity of c , and combining this with (|4.29|) . It follows from (|4.33p that the term 
J r OZz-vdx in the momentum equation (|4.3ip . with a selection z G L°°(r c x (0,r)) from d(dtu) as 
in "the second of ([4~3"T) . is well-defined for every v G W. Observe that the L°°(0, r; L 2+V (T C ; ir- 
regularity of \i derives from (|3.13p . 

To prove uniqueness, we proceed as in [Til Sec. 5]: given two solution pairs (ui, /u x ) and 
(u2, fi 2 ) to the Cauchy problem for (|4.3ip . we test the equation fulfilled by u := Ui — u 2 and 
fi := fi x — fi 2 , Hi = |K(</>g(uN)n)|zj for some z^ G d(d t u.;), i = 1,2, with v := d t u and integrate 
in time. With straightforward calculations, setting z = zi — z 2 and using the place-holder 01 from 
(|4.33p . we obtain 



b(d t u,d t u) da + -o(u(t),u(t)) 
o 2 



< 



Xu • d t u dx ds 



o Jr c 



Hz ■ dtudxds 



6^((ui) N )-^((u 2 )N))n • dtudxds 



Now, the third term on the left-hand side is positive by monotonicity of d and positivity of c (cf. 
(|3.9p ). whereas the second integral on the right-hand side can be estimated relying on the Lipschitz 
continuity of <f>' e . All in all, the desired contraction estimate for u ensues from applying the Gronwall 
lemma, with calculations analogous to those in [TTJ Sec. 5]. Clearly, from this we also deduce that 
there exists a unique \x complying with (|4.31[) . 
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Finally, to prove estimate (|4.32p we test (|4.3ip by v = 9 t u and integrate on (0, t) with 
t € (0, r] . We obtain the following estimate 



C h [ ||9 tU ||^ds+^||u(i)||^+ f f c(d-d s )n-d t udxds+ I <j> £ (Mt))dx 

< c^||uoIlw + ^ <t>e(MQ))te + jf (0\\h + ||F|| W ')IMwds + ^ ||X|| ffrc ||u||w|ia t u|| w d.s 

< c(l + M 2 + £ ||X|| 2 ffrc \\u\\ 2 w ds) + yII 9 * u IIw, 



where for the left-hand side we have used the chain-rule identities 
ft 



r 11 

J a(u(s),d t u(s)) ds = -a(u(t), u(t)) - -a(u , u ) 

/ / <^(u N )n ■ dtudxds = / (/) E (u N (t)) dx - 
Jo Jr c Jr c 



(4.34) 



(4.35) 



<j) £ (u-u{0))dx, 



(cf. [TBI Lemma 3.3], [18j Lemma 4.1]) the coercivity properties (|3.T[) and (|3.8[) of the forms a and 
6, and ultimately that J T </> e (mn(0)) dx < y(u ) < oo by (|3.23p . By the monotonicity of d and the 
positivity of c (cf. (|3.9p ). we also infer that 

/ j c(d-d s )(j, -d t udxds > 0. (4.36) 
Jo Jr c 

For the estimates on the right-hand side of (|4.34|) . we rely on Young's inequality and trace theorems. 
Therefore, (|4.32[) ensues from (|4.34[) by the Gronwall lemma. Clearly, the estimate for z follows from 
the fact that |z| < 1 a.e. on T c x (0,r), whence the estimate for fi , in view of (|3.13[) . ■ 



Thanks to Lemma 14.61 we may define an operator 

7 X :M T ^U T ^{ueif^O^W) : ||u|| H i ( o iT; w) < (4.37) 

mapping every quadruple (i?,i? s ,u, X) £ y T into the unique solution u (together with fi e |3l(0 E (wN)n)|d(9tu) ) 
of the Cauchy problem for (|4.31[) . 

Step 2: As a second step, we solve (the Cauchy problem for) (|3.44p - (|3.45l) . with $ s £ tt2(^t) and 
u from Lemma 14.61 on the right-hand side of (13.441) . 

Lemma 4.7 Assume (|3.ip . Hypotheses (IV)-(V), and suppose that Xq complies with (|3.24l) . 

Then, there exists a constant S2 > such that for all ($ s ,u) £ ^2(^1-) x 11 t there exists a 
unique pair (X,£) £ (L 2 (0, r; H 2 (T C )) n L°°(0, r; Vr c ) PI i? 1 (0, r; i?r c )) xI 2 (0,T;iJ r J fulfilling the 
initial condition (|3 .38|) . </ie relations 

d t X + AX + i + a'(X) = -\'(X)d s -^\u\ 2 a.e. mY c x (0,T), 
£ £ /3(X) a.e. m T c x (0,T), 

and £/ie estimate 



II^IU 2 (0,r;H 2 (r c ))nL=-(0,r;y rc )nH 1 (0,-r;Hr c ) + II £ II L 2 {0,t;H Tc ) ^ ^2- (4.39) 

Proof. The well-poscdness of the Cauchy problem for (|4.38|) follows from standard results in the 
theory of parabolic equations with maximal monotone operators, after observing that, by Sobolev 
embeddings and trace theorems, & s and 1/2 1 u| 2 arc respectively estimated in L 2 (0, r; L 2 / S (Y C )) 
(with 5 £ (0, 1) as in (|4~2^|) ). and in L 2 (0, r; H Tc ) ■ 
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In order to obtain estimate (|4.39p . we test (|4.38|) by dfX and integrate in time. Exploiting 
the chain rule for /? from [TBI Lemma 3.3] 



we obtain the estimate 
rt 



o Jr 



1 



£d t Xdxds = I P(X(i))dx-J P(X )dx, 



(4.40) 



o Jr c 
1 



\d t X\ 2 dxds + ^ j \VX(t)\ 2 dx + I '0(X(t))dx 



< 



|VX | dx 



P(X )dx 



(4.41) 



<r'(X)d t Xdx ds 
<C{\\X \\ 2 Vr +l) + h+h + h 



JT € 
2 



\'(X)d s d t Xdxds 



o Jr 



\u\ 2 d t Xdx ds 



o Jr c 



Now, the last inequality ensues from (|3.24[) , and we estimate the integral terms , i = 1,2,3, as 
follows: 



\d t X\ 2 dxds + C 



o Jr c 
t 



(|X| 2 + l)da;ds + C" 



< 



4 Jo 



JT C 
t 



d t x\ 2 dxds + c / \\d t x\\l 2[0jS . Hrc) ds + c\\x \\ 2 Hrc + c' 



h<c {\\X\\ L 2ni-s) { r c ) + l)\\tis\\LVHr c )\\dtX\\Hr c ds 



(4.42) 



h < 



<7 \\dtX\\ 2 Hrc ds + C \\Mi»/i(r )(\\X\\vr c + 1 ) ds > 



\\d t Xf H ^ds + c I \\u\\ 4 w ds, 



\H Tc 



where for Ii we have used that a' is Lipschitz continuous by virtue of (|3.16[) . whereas the estimate for 
I2 follows from the Lipschitz continuity of A' and from observing that i/ 1 (L c ) embeds continuously 
in i 2 /( 1_<5 ) (r c ) for any S G (0, 1), and for ^3 we have exploited the trace result (|3.4[) . As for the 
left-hand side of (|4.4ip . it remains to observe that, by convexity, 



p(X(t))dx>-c' \X{t)\dx-C. (4.43) 



3 c', C > such that for every t e (0, r] 



Plugging (|4.42l) and (|4.43[) into ()4.41[) and using the Gronwall Lemma we conclude an estimate for X 
in L°°(0, r; Vr c ) H H 1 (0, r; i?r c ) ■ A comparison argument in (|4.38[) (cf. also the forthcoming Seventh 
a priori estimate in Sec. I4.3p . then yields an estimate in L 2 (Q, r; -Hr c ) for AX + £ , hence for AX and 
£ separately in L 2 (0,t; flr c ) by monotonicity of /3. Therefore, by elliptic regularity results we get 
the desired bound for X in L 2 (0, t; H 2 (T C )) , which concludes the proof of (|4. 39[) . ■ 

It follows from Lemma 14.71 that we may define an operator 



J 2 :i 2 (WxU r ^l r :={XeL 2 (Q,T;H 2 (r c ))nL ac (0,T;H 1 (r c ))nH 1 (0,T;H rc )nX T 

IMU 2 (U-r;ff 2 (r c )) + |M|i°°(o,T;VrJ + IMIffMtVnflrJ ^ ^2} 



(4.44) 



mapping (i? s , u) G ^Cdr) x U t into the unique solution X of the Cauchy problem for (|4.38[) (together 
with £ G /3(X)). 



Step 3: Eventually, we solve the Cauchy problem for the system (|4.25[ l4~26|) with hxed ($,i? s ) G 
Ti,2 (Vt) and (u, X) from Lemmas 14.61 and 14.71 respectively. In particular, we set 



J := fc(X)(tf - 8 ) + c'(tf - ^)|^(^( UN )n)||a t u T | 



(4.45) 
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and plug it into the boundary integral on the left-hand side of (|4.25[) and on the right-hand side of 
(|4.26|) , Observe that, due to (|4.30|) . (|4.32|1 . (|4.39|1 . and to the Lipschitz continuity of c and k, there 
holds 

S G L 2 (0, r; L 4/3+s (r c )) for some s = s(S) > 0. (4.46) 

We mention in advance that, relying the very fact that the boundary term 3 in (|4.47p below does not 
depend on the unknown $ , we will be able to prove uniqueness of solutions for (the Cauchy problem 
for) dQ7t . 

Lemma 4.8 Assume ()3. 1|) . Hypotheses (I), (II), (III), and (V), suppose that h complies with (|3.17|) . 
and let (i?§,^ e ) fulfill (14441) . 

Then, there exist positive constants S3 and 54(e) such that for all (i9,i? s ,u,X) G Trinity t) x 
tl r x X r i/iere exists a unique couple of functions i? s ) complying with (|4.19p - (|4.22[) . fulfilling the 
initial conditions (|4.23[1(|4.24|) . £/ie equations a.e. in (0,T) 

(d t £, e ($),v) v - [ div(<9 t u)uda; + / WVudx + [ 3vdx = (h,v) v VtigF, (4.47) 
Jn is! Jr c 



{d t L e {d s ),v) Vr ~ d t X(X)vdx + S7$ s Vvdx = Jvdx Vu G Vr c , (4.48) 
and i/ie estimates 

H^llL 2 (0,r;V)nL~(0,T;Li(O)) + II II L 2 (0,r ;V Ta )nL°° (0,r ji 1 (r c )) < S3, (4.49) 

||9 t £ e (i9)|| L2(0:r;yo + \\d t Z e (# s )\\ L z {0tT . v , a) < S 3 , (4.50) 

ll*IU-(0,r;fl) + ||0.||i-(O,Ti!rr o ) < (4.51) 

Proof. Observe that system (|4.471 I4.48[) is decoupled, hence we will tackle equations (|4.47[) and 
(|4.48p separately. 

Also taking into account (|4.46|) . the wcll-posedness for the Cauchy problem for the doubly 
nonlinear equation (|4.47[) follows from standard results, cf. QI51 Thm. 1] (see also [§J Lemma 3.5]): 
in particular, uniqueness for (|4.47p is trivial, since the terms div(<9 t u) and 3 are fixed. In order to 
conclude estimates l[Og]l - ([4T5T]l for , we test (|Q7|) by -d and integrate on (0, t) with t G (0, r] . 
Recalling the definition (|4.6[) of 3 e , we exploit the formal identity (cf. Remark 14.121) 

(3^(0), 0> v = / 4(tf)d t iMdai = 4: f JeWda (4.52) 
Jn d* Jn 

and thus infer 



| / |??(i)| 2 dx + Ci / |0(t)|dx-C 2 + / / |Vi9| 2 dxds 
2 Jn Jn Jo Jn 



< f j £ ($(t))dx+ f f |Vtf| 2 d.Tds 
Jn Jo Jn 

< f j £ (# s )dx+ f f div(5tu)(0-m(0))da;ds- / / - m(0)) dxds 
Jn Jo Jn Jo Jr c (4.53) 

+ / (h,d-m(d)) v ds+ [ [ div(d 4 u)m(tf) dxds 
Jo Jo Jn 

3m('d)dxds+ f f hm(§)dxds 
o Jr c Jo Jn 

< S^l + ||0 o |Ui(n)) + h+h+h + h + h + h, 

where the first inequality is due to (|4.8b|) . and the estimate for J n J e (^o) ^x follows from (|4. 14e[) . 
As for the terms li , i = 1,...,6, by the Sobolev embeddings and trace results (|3.4j) . joint with 
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Poincare's inequality p.5[) , we easily have 



ds 



h+h+h< j (||Mlw||l?-m(i?)||H + ||^| i 4/3 (rc) || 1 ?-m(7?)|| L . ( r c ) + ||/i||vHI^-^WIIv 
< \J o ^ + j (\\d t u\\ 2 w +\\ n 2 L * /HTc) + \\h\\v>) ds, 

h + h+h<C j (||MIw+||^/3(r )+IWU»(n)) ¥\W{U) ds. 

We plug the above estimates into (|4.53[) and use (|3.17[) . estimate (|4.32p , and (|4.46[) . Relying on the 
Gronwall lemma, we conclude estimates (|4.49l) and (|4.51[) for -d . Estimate (|4.50|) for £ s (#) follows 
from a comparison in ()4.47j) . 

Since the calculations related to the analysis of equation (|4.48|) are completely analogous, we 
choose to omit them. ■ 



Thanks to Lemma 14.81 we may define an operator 
T3 : 7ri l2 (V T ) x U T x X T -> 

W T :={(*, 1?.) 6(L 2 (0,T;y)ni m (0,T;fi)) x (L 2 (0,t; VrJ n L°°(0,r; ifrj) : 

(4.54) 

II^IU 2 (0,T;V)nL^(0,r;L 1 (n)) + 1 1 $s \\ L 2 (0,r ; Vr c )ni°° (O.t;!, 1 (T c )) < ^3, 



= (0,r;H) + \\^s\\L^{0.T-H rc ) < S^e)} 



mapping (t?, $ s ,u, X) € tti^C^t) x U r x X r into the unique solution of the Cauchy problem 

for system (j4~47l |448)) . 

We are now in the position to prove the existence of local- in-time solutions to Problem (P e ) , 
defined on some interval [0, T] with < T < T. Note that T in fact does not depend on the 
parameter e > . 

Proposition 4.9 (Local existence for Problem (P £ )) Assume (|3.ip . Hypotheses (I)-(V), and 
conditions 1(347]) - ([349]) on the data h, f, g, ([333]) - pTSI]) on u , X , and 1444]) on tf°< £ . 

Then, there exists T 6 (0,T1 suc/i £/ia£ /or every e > Problem (P s ) admits a solution 
i? s , u, X, /x , £) on the interval (0, T). 

Proof. Let the operator T : y T — > W T x U r x I T be defined by 

T(tf,£ s ,u,X) := (tf,i? s ,u,X) with J X := T 2 (£ s , u), (4.55) 

[(0,0.) :=T 3 (^^,u,X). 

In what follows, we will show that there exists T G (0, T] such that for every e > 

T maps into itself, (4.56) 
T : — > is compact and continuous w.r.t. the topology of 
L 2 (0,T;fr 1_ '(n)) x L a (0,T;H 1_a (T c )) x L 2 (0, t; i/ 1_<5 (il; M 3 )) x L 2 (0,T;H r J. 



(4.57) 



Ad (|4~56]) . In order to show ([436]) . let (i?,i? g) u,X) 6 ^ T be fixed, and let (i?,i? s ,u,X) := T(i?,i? s ,uX). 
We use the interpolation inequality 

l|^WII^-'(n) < c|W*)lljfi(n)IW*)lli'(n) for a.a.t e (0,r) (4.58) 

(cf. e.g. [TTl Cor. 3.2]). Now, a further interpolation between the spaces L 2 (0,t;V) and L°°(0, t; L 1 ^)) 
and estimate (|4.49|) also yield the bound ||i?||iio/3(o t-l 2 (si)) ^ &S3 for some interpolation constant 
C . Integrating (|4.58p in time and using Holder's inequality we therefore have 



\ 2 L W i-* m <cj o II^WII^ 11^)11^(0)^ (4 59) 

- C \\ V \\l 2 (0,t-H 1 (Q,)) Z ll W llL 10 /3(0,T;L 2 (f2)) - ^^l 1 
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We use (|4.58p for $ s and u to perform calculations analogous to (|4.59p , whereas for X we trivially 
have ||X]||,2(o ph t ) — ^IMIi°°(o t-H r ) — ^2- Combining all of these estimates, we conclude that 
there exists a sufficiently small T > for which (|4.56|) holds. 

Ad (|4.57j) : compactness. Exploiting estimates (|4.32[) . (|4.39[) . and the compactness results [30l 
Thm. 4, Cor. 5], it is immediate to check compactness of the operator T as far as the (u, X) -component 
is concerned. As for the (#, i? s ) -component, from (|4.49[) - (|4.5ip and the Lipschitz continuity of L e 
we deduce an estimate (with a constant depending on e) for £ e (i9) in L 2 (Q,T;V) n -L°°(0, T; if) n 
^{O^V) (for L e (d s ) in L 2 (0, T; VbJnL°°(0, T; HrJnH 1 ^, f; V/J , resp.), whence compactness 
for £ e (i?) in L 2 (0,f; H 1 ' 8 (Q)) (for £ e (tf s ) in L 2 (0, f ; i/ 1 -" 5 ^)) , resp.), hence for i9 = ^(^(tf)) 
in the same space (and analogously for i? s ) . Observe that the latter argument relies on the bi-Lipschitz 
continuity (|4.9[) of L e . 

Ad (|4.57[) : continuity. In order to prove that T (|4.55|) is continuous, we will check that the operators 
Tj, i = 1,2, 3 defined by (|4.37p . (|4.44p . and (|4.54p are continuous w.r.t. to suitable topologies. 

First of all, we fix a sequence {(t? n , $ s>n , u„, X„)}„ C converging to a (i9oo, ^ S;00 , Soo, Xoo) £ 
, with 

ini 2 (0,f;ff u (O)), ?a,n->?.,oo inX^O.fjfr 1 -*^)), 

u„^u oc inL^O,?;^ 1 - 5 ^;^ 3 )), X n -> X^ in L 2 (0, f ; ff r J 

as n — > oo. We let u„ := Ti($„, t9 Si „, u„, X„) , and denote by (/jj„ the associated sequence such 
that (u„, fi n ) fulfill (|4.3ip . Due to estimate (|4.32p . there exist a (not relabeled) subsequence and a 
pair (UqojZoo) such that as n — > oo 

u„^ Uoo inlT 1 (0 j r;W) 1 Mn^>oc in £°°(0, T; L 2 +"(r c ; M 3 )) . (4.61) 

Hence, by well-known compactness results, (u„)„ strongly converges to u in C°([0, T]; i7 1_5 (f2; R 3 )) 
for all 5 G (0, 1] . Now, combining convergences (|4.60p and (|4.6ip and arguing in the very same way 
as in the proof of Thm.Q](cf. the forthcoming Scction[5]), we manage to pass to the limit as n — >• oo in 
(|4.3ip . concluding that the pair (u^, fi^) fulfill equation (|4.3ip with (??oo, $ s ,oo, UooXoo) . Therefore, 
we have that 

Uoo = Ti^ooj^oojUoojXoo), and convergences (|4.6ip hold for the whole {(u„, £t„)}„, (4.62) 

the latter fact by uniqueness of the limit. 

Secondly, we consider the sequence X„ := ^{fls.n, u n) with (u„)„ from the previous step, 
and let (£ n )n be the associated sequence of selections in j3(X n ), such that (X„,£„) fulfill (14.381) . 
Thanks to estimate P~39"]) . we have that (X n ,£„) n is bounded in (L 2 (0, f ; H 2 {T C )) n L°°(0, f ; VrJ n 
i? 1 (0, T; Hr c )) X L 2 (0, T; Hr ) ■ Therefore, there exists (Xoo, £oo) such that, up to a subsequence, as 
n — > oo 

X„-*X QO in L 2 (0, T ; ii 2 (r c )) n L°°(0, f ; V Te ) n ^(0, T ; iT r J, 



^^Coo inL 2 (0,T;ff : 



(4.63) 



and (X„)„ strongly converges to Xoo in L 2 (Q, T; ii 2 ~ p (r c )) n C°([0, T]; H 1 ^ <5 (r c )) for all p G (0,2] 
and <5 G (0,1] by [30j Thm. 4, Cor. 5]. Relying on convergence (|460|) for ($ s , n )„, (|4~6Tj) for (u„)„, 
and (|4.63p , and arguing as in the passage to the limit developed in Sec. 03 it can be shown that the 
functions (X,£) solve (|4.38[) with (^ s ,oo, u oo) , i-e- 

Xoo = 72(i? Sj c») u oo)i and convergences (|4.63|) hold for the whole (X n , £ n )„ . (4.64) 

Thirdly, we let (i?„,i? s ,„) := T 3 ($„, i? S) „, u„, X„) with (u„)„ and (X„)„ from the previous 
steps. Estimates (|4.49p - (|4.5ip imply that there exist ($oo, $s,oo) such that, along a (not relabeled) 
subsequence, the following weak convergences hold as rt — » oo 



tfn-^oo mL 2 (0,T;V)nL°°(0,T;H), 
^ s ,n^*^ l00 in L 2 (0, f ; K r J n i°°(0, f ; H Te ) 



(4.65) 
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Furthermore, taking into account the that (<C £ (i? n )) n ((£>e(^s,n))n > respectively), is bounded in 
L 2 (0 > f;V)nL°°(0,f;H)nH 1 (0,f; V) (in L 2 (0, f ; V Ta ) nL°°(0, f ; H r J nH\0, f ; V^J , resp.) and 
relying on [501 Thm. 4, Cor. 5], we find that £ e (0 n ) -> £ e (i?) in L 2 (0, f ; H 1 - 5 ^)) n C°([0, f ]; # ) , 
and analogously for (£ e (i? s ,„))n • Therefore, thanks to the bi-Lipschitz continuity of L e we conclude 
that 

n -H?oo m£ 2 (0,T; 1^(0)) nC°([0, ?];#), 
^ s ,n -> in ^ 2 (0, f ; ff 1-4 (r c )) n C°([0, ?]; HrJ. 

Convergences (|4T60|) for ($ n , d s . n ) n , (|4~6T]) for (u„)„, (|4~63)l for (X„)„, and (|4~65|) - (l466l) . combined 
with the arguments of SecJH allow us to pass to the limit as n — > oo in system (|4.47[ 14.451) . Therefore, 
we conclude that 

(tfoo,$ s ,oo) = ^(^ooj^oojUoojXoo) and 

(4.67) 

convergences (|4.65p - (|4.66l) hold for the whole (^n, $s,n)n- 
Ultimately, the continuity of T ensues from (|4.62p . (|4.64p . and (|4.67l) . ■ 

Remark 4.10 A key trick to prove compactness and continuity of the operator T (|4.55[) in the 
(■$, ■$<,) -component has been: 

- to prove compactness and continuity in (£ £ ($), L E (# s )) (exploiting the estimates on the pair 
(<9 t £ e (i9), dfCei'&s)) deduced from a comparison in the temperature equations (|4.25|) and (|4.26p ). 

- then to infer compactness and continuity in # s ) , relying on the fact that L £ is bi-Lipschitz, 
cf. Lemma l4~2l 

Obviously we will not be in the position to use such an argument any longer, when taking 
the limit as e — >• 0. Indeed, for such a passage to the limit we will rely on new, BV-type estimates 
on # s ) (cf. the Fifth and Sixth a priori estimate in the forthcoming Sec. I4.3p . Exploiting such 
bounds and a version of the Aubin-Lions theorem for the case of time derivatives as measures (see, 
e.g., Chap. 7, Cor. 7.9]), we will conclude the desired compactness for •& and # s . 

4.3 The approximate problem: global existence 

We now prove the following 

Theorem 4.1 (Global existence for Problem (P e )) Assume (|3.ip . Hypotheses (I)-(V), condi- 
tions (|3.17p - (|3.19p on the data h , f , g , (|3.23p - Q3.24p on Uo , and Xo , and (|4.14p on the approximate 
data $o an d ^s' e • Then, 

1. for all e > Problem (P £ ) admits a solution i? s ,u, X, on the whole interval (0,T); 

2. there exists a constant C > such that for every e > and for any global-in-time solution 



(#, i? s ,u, X, /i,£) to Problem (P e ) the following estimates hold: 

£ 1/2 ||tf||L~(0,T;H) +£ 1/2 ||l?«||L»(0,r ; ffr c ) < C > ( 4 - 68 ) 

ll#l|L 2 (o,T ; v)n£ c >°(o,T;z, 1 (n)) + ll^s||L 2 (o,T;y rc )nL°°(o,T;L 1 (r c )) < C s (4.69) 

ll u llii' 1 (0,T;W) + ||^]|l,»(0,T;Vr c )ni? 1 (0.' r ;' H 'r c ) - (4.70) 

ll'? !> £( u N)n|| L 2 (0iT; Y^) + ||/*||L«(o,T;L*+"(r c ;R!»)) < C wrf/i v > /rom P-13P, (4.71) 

||a^C e (i?)IU»(o,r i V') + l|9 t £ E ^ s )||L 2 (o,T;Hi(r c )') < C, (4.72) 

II^I|bv(o,t ; wi.9(o)') < C for any q > 3, (4.73) 

ll^s||BV(o,T;W 1 .>'(r c )') < C / or any a > 2, (4.74) 

ll X IU 2 (0,T;_ff 2 (r c )) + ll^llL 2 (0,T;_ffr c ) ^ ( 4 - 75 ) 

l|£ 6 Wlli-(o,Tifl) + ll^s(^)IU-(o,T;H rc) < C . (4.76) 



In order to prove Thm. 14. li in what follows we establish a priori estimates on the (i?, i? g ,uX)- 
componcnt of any given solution to Problem (P e ), independent of the time-interval on which such 
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solution is defined. Exploiting these global-in-time estimates and a standard prolongation argument, 
we will conclude that the local solution to Problem (P e ) from Proposition ^. 91 extends to a global-in- 
time solution. In this way we will obtain the first part of the statement. 

In fact, it will be clear from the calculations below that such global estimates hold for a 
constant independent of the parameter e > 0, whence (|4.68p - (|4.76p . which will provide the starting 
point for the passage to the limit as e — > in Sec. [5J 

We mention in advance, the First, Fifth, Sixth a priori estimates below are only formally 
derived: in Remark 14 . 1 21 later on we will clarify how they can be made fully rigorous. 

Notation 4.11 We stress that, from now on the symbols c, c', C, C , shall denote a generic constant 
possibly depending on the problem data but not on e, which we let vary, say, in (0, 1). Since the 
estimates below are not going to depend on the final time, we will perform all the related calculations 
directly on the interval (0, T) . 

First a priori estimate. We test gjggl by 0, flOg) by S , P~2"7|) by d t u, and ([PU) by d t X, 
add the resulting relations, and integrate on (0,t), t £ (0, T] . Recalling (|4.52[) . we formally have 

' {d t L £ (#),#} v ds= f J e {-d{t))dx- [ 3e($o~)dx 
o Jn Jn 

> 1\W)\\h + CiWmhnn) - + \\*oh H a)) ( 4J7 ) 

>l\Wt)\\% + Ci\W)hw-c, 

the first inequality due to (|4.8bj) and (|4.14e|) . and the last one to the first of (|3.21j) . In the same way, 
we have t 

(dt& B (4.Wvr e ds ^ fll^Wllk +Ci\\Mt)\\Li(r o) -C. (4.78) 

We take into account the cancellation of some terms, the chain-rule identities (|4.35[) and (|4.40[) . as 
well as 

j J (xu-d t u+^\u\ 2 d t X^ dxdr^^J X(t)\u(t)\ 2 dx - ^ J X |u | 2 da;. 
Therefore, with easy calculations we arrive at 

£ ""'-"2 , n ILa/Vlll ....if f Y7„Q|2 i j„ , e i|„Q l.\\\2 



\\m\\B+ci\\mh^) + / / m'dxdr+-\\ut)\\H rc +ciiitf s (*)iu 1( r c) 



\Vd s \ 2 dxdr + f I k(X)(tf-tf s ) 2 dxdr 
Jo Jr c Jo Jr c 

+ [ f c'(ti - $ S )(V - i9 s )m'Au N )n)\\d t u T \dxdr + [ b(d t u, d t u) dr 
Jo Jr c Jo 

+ io(u(t),u(t)) + i / X(t)\u{t)\ 2 dx + f f c(ti-tf s )n-d t udxdr 
+ / (f> e (u N (t))dx + 



(4.79) 



I I \d t X\ 2 dxdr + \ [ \VX(t)\ 2 dx 
Jo Jr c 2 J Tc 

< C + ia(u ,u ) + ^ X |u | 2 d.T + J e («N(O))da;+i||VXo|||- r(> +h + h+h + h- 

Now, observe that, due to the positivity of k in fl3.14p . the seventh integral term on the left-hand side 
is positive, and so are the eighth term, thanks to Q3.9p . the eleventh, since X £ dom(/3) C [0, +00) , 
and the twelfth, by (|3.9j) and the fact that fi ■ dtu > a.c. in T c x (0, T) . As for the right-hand side 
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2 

L 2 (0,T;W) 



of (|4.79[) , it holds J r </> e (un(0)) dx < <p( u o) < oo by (|3.23[) . Moreover, we have that 
h=J (F,d t u) w dr<g^ ||9tu||^dr + C||F||2 2(0T . w0 < ± b(d t u,d t u) dr + C\\F\ 

h = I {h,§) v dr= [ {h,ti-m{-&)) v dr+ [ [ hm{d)dxdr 
Jo Jo Jo Jn 

<\ ! [\W\ 2 dxdr + c[ ||&||jr||0|| L i C n)dr + c/ \\h\\ 2 v ,dr, 
1 Jo Jn Jo Jo 

where the second inequality in the first line is due to (|3.8j) . and we choose g = Cb/2 in order to 

absorb the term J* b(dt\i, dtu) dr into the corresponding term on the left-hand side. In the estimate 

for I2 , the second passage follows from Poincare's inequality. Finally, exploiting the convexity of (3 

and the fact that a' has at most a linear growth, (cf. also (|4.42l) and (|4.43[1 ). we have 

h = ~f P(X(t))dx<C f \X{t)\dx + C <- f ( \d t X\ 2 dxd S + C\\X \\ Llin) +C>. 



h = - J jf a'{X)d t Xdxds <lj o J r \d t X\ 2 dxds + cj* |P t X||| 2(0 ^ rJ ds + C'\\X \\ 2 Hr 



C. 



We plug the above calculations into the right-hand side of (|4.79[) . Relying on assumptions (|3.17[) for 
h, (|3.20p for F, and on (|3.23p - (|3.24[) for the data Uo and Xo, applying the Gronwall Lemma we 
immediately deduce estimates (|4.68|l - (|4.70|) . 



Second a priori estimate. It follows from estimate (|4.70[) and the continuous embeddings ()3.4|) 
that the term Xu on the left-hand side of P~2"7)l is bounded in L 2 (0, T; i 4_e (r c ; R 3 )) for every 
e £ (0, 3] . Therefore, taking into account the previously obtained estimates on u and 1} , and arguing 
by comparison (|4.27|) . we also obtain the bound 

||C(0 + ^(«N)n|| L2( o,T;Y^) < C. (4.80) 

Exploiting the fact that \x and (j) f e (u^)n are orthogonal and arguing as in Sec. 4], from (|4.80|) 
we conclude that 

||c(l?-^)/i|U 2( o,T ; Y^) + ||^("N)n|U 2(0 ,T;Y^) < C, (4.81) 

whence (|4.71| : the estimate for /x = |3i(</4(uN)n)|z follows from the fact that 01 : L 2 (0,T; Y' r J -t 
L°°(0,r;L 2+,y (r c ;]R 3 )) is bounded thanks to ([57lt| . and from the fact that 

|z| < 1 a.e. on T c x (0,T) (4.82) 

by the definition Q1.13P of d . 



Third a priori estimate. We argue by comparison in the temperature equation (|4.25|) . It follows 
from estimates (|4.69p - (|4.70[) . the continuous embeddings (|3.4[) . and the Lipschitz continuity (|3. 14[) of 
k, that 

||fc(X)(7?-i9 s )|| L2(0iTiL1 _« (rc)) <C for every ee (0,3]. (4.83) 

Now, taking into account estimates (|4.70|) and (|4.71[) . again (|3.4[t . the fact that 01 : L 2 (0,T; Y' r ) -> 
L°° (0,T; L 2+v (r c ;R 3 )) is bounded by (|343]) . and the boundedness of c' by (|3^|) . we have at least 

||c'(i? - ^)|^(^( UN )n)p t u T ||| i2(0iT;L4/3(rc)) < C. (4.84) 

Therefore, in view of ()3.17j) on h, by comparison in (|4.25| we obtain estimate (14.72j) for <9i£ e ($). 

Fourth a priori estimate. We rely on (|4.83| . (|4. 84[) . and estimate (|4.70j) which, combined with 
(|3.15[) for A , in particular yields 

\\dtX(X)\\ L 2 (0iT . L3/2(To)) . (4.85) 
Therefore, a comparison in the temperature equation (|4.26|) yields the second of (|4.72|) . 
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Fifth a priori estimate. Let us formally rewrite (|3.39p as 
£' e (d)d t tf ■ v dx = [ div(<9 t u) vd.T - / VtfVvdx- [ fc(X)(t? - -& s )v dx 



c'($-'3 s )\3l((j)' e (u N )n)\\dtU T \vdx+ {h,v) v Vv G V a.e. in (0,T) 



(4.86) 



and choose in (|4.86p a test function v G V of the form 

1 



with w G and g > 3. 



(4.87) 



Taking the contraction property f|4. 10[) of into account and considering that ?9 G V C L 6 (5!) , we 
have that jTjjpj G V and therefore £ , 1 ^. ) u> G . Furthermore, it follows from (|4.10[) that 



< 



1 + ep e (l) + e 2 



I'd- 1| < |t?| +2 a.e. inQ x (0,T), (4.88) 



where we have also used formula (|4. 13[) involving the resolvent p e of In, which satisfies p e (l) = 1. 
Therefore again exploiting (|4. lOf) we find 



^(tf) 



< l|V^||ff 



Now, we have 



1 



-to ) dx = / dt"dwdx. 



(4.89) 



(4.90) 



Moreover, in view of (|4.10[) , (|4.88|) . (|4.89p . the previously obtained estimates, as well as (|3.17[) on h, 
we see that 



dx 



< liatullwlKlltfHfl + c)HU~(n) = /i ei 2 (0,T), 
< l|Vi>||i||t«||ioo ( n) + (||<?|U» ( n) + c)||Vi?||jr||Vu;|U« ( f l ) =/j£ ^(O.T), 



< 



p(X)(0 - 0.)llffr e (l|0|U«(r e ) + c)||io|U*CTo) = /3 G ^ 2 (0,T), 



< 



c'(tf - ^)|^(w N )n)||5 t u T |— — -wdx 
||c'(tf - ^)|K(^( U N)n)||a t UT||| L 4/3 ( r c )(||^IU4(r B ) + l)IMU~(r.) = h e ^(0,T), 



w > 



< 



|v ((IMIv + c)|M|i- ( n) + (||0||i«(n) + c)||V«;|| L 3 (n) ) = f 5 G ^(O.T), 



where we have also used the continuous embeddings (|3.4p . W 1,q (il) C L°°(f2), V C i 6 (51), as well 
as the trace result W 1,q {£l) C L°°(r c ) . All in all, we conclude that 



Hence, we have 
yielding (|4T75|) . 



i/ei^T) for a.a. i € (0, T) V«J G W 1,Q, (f2) : 

d t ${t)wdx = (d t 'd(t),w) wl , qin} < f{t)\\w\\ w i, q{n) . 



(4.91) 



(4.92) 
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Sixth a priori estimate. Wc proceed in an analogous way with p.40[) and test it by 

-w, with w G W 1,a (T c ) and a > 2. (4.93) 



The analogues of estimates ([4.88)1 and ([4.89|) hold, therefore we obtain 



= h+h + h + h- 

(4.94) 

Using that i and A' are Lipschitz and relying on the continuous embedding M /1,CT (r c ) C L°°(r c ) 
we have 

|4| < ||A'(X)|U4 (rc) ||S t X|| L2(rc) (||^|| i4(ro) + c)|MU-<r ) 

< C(l + !|X|| L 4 (ro) )||a t X|| i2(rc) (||^|| i 4 (rc) + c)||HU~(r o) = h e L\0,T) 
in view of estimates ([4.69)1 — (|4.70j) . Analogously, we have 

\h\ < ||W s ||ffr c |M|i~(r ) + (ll*«IU'(r.) + c)\\V# s \\Hr c ||V«;||i-(r ) = fr e ^(O.T) 

IM < (l + IIXlKJII^-^IU^dl^lU^+c)!^!!^^) i/sSl^T), 

|4| < ||c'( l ?-^)|^(uN)n)||0tUT||| i v3 ( r c) (||^||L4 ( r c )+c)|| W |Uoo (rc) = / 9 eL^O.T) 

where in the estimate of Jg we choose p in such a way that 1/p + 1/2 + 1/er = 1, exploiting (|3.4p , 
whereas to deal with 4 we have used the fact that k is Lipschitz, and finally for Ig we have relied 
on ([4.84)1 . All in all, we conclude that 

3/etf(0,T) for a.a. t G (0, T) Viu G W 1,ff (r c ) : 

(4.95) 

<9 t tf s (i)u;dx = (dttf s (t),w} w i, a(r ) < f{t)\\w\\w 1 '<'(r c )- 

r c 

Hence, we have 

\\d t # s \\Li(o,T;W^(r y)<C, (4.96) 

whence (j4.74|) . 

Seventh a priori estimate. We test (13.44)) by AX + £ and integrate in time. Taking into account 
the chain-rule identity ()4.40p , we obtain 

i / \VX(t)\ 2 dx+ f /3(X(t))dx + f f \AX + £,\ 2 dxdr 



r c Jr c Jo Jr c 

|VX | 2 d2;+ / p(X )dx + I 9 + Iio + 111, 



1 

~~ 2 
and estimate 

h = - I I a'{X)(AX + £,)dxdr<c[(l + \\X\\ Hrc )\\AX + Z\\ Hrc dr 
Jo Jr c Jo 

<\J \\AX + i\\ 2 Hvc dr + cj \\X\\ 2 Hrc dr + C, 

ho = - f f \'{X)d s {AX + t)dxdr<C f {l+\\X\\ LP{Ttt) )\\V s \\ LnTc) \\AX + (i\\ Hrc dr 
Jo Jr c Jo 

1 '* 



'ii 



<C(l+||X|Uo. (0 , Ti vi. o )ry ¥s\\v Tc dr +l I \\AX + a 2 Hrc dr 
\f j \u\ 2 {AX + t)dxdr< 1 - f \\AX + i\\ 2 Hr€ dr + C f M^dr, 
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where for Ig we have used the Lipschitz continuity of a' , for iio chosen p and v in such a way 
that 1/p + l/v + 1/2=1 and then used the continuous embedding (|3.4j) . and analogously for Jn . 
Applying the Gronwall Lemma and taking into account estimates (|4.69p - (|4.70p . we then conclude 
that J Q \\AX + £||jy r dr < C. A monotonicity argument yields 

||AX|| i2(0 , T . HrJ + ||e|| ia(0l T i ffr <> ) < C > ( 4 - 97 ) 

whence (|4.75p by elliptic regularity. 

Eighth estimate. We test (|4~25|) by £ E (i?) G V and ([4~26)l by £ e (tf s ) G V Tc , add the resulting 
relations and integrate in time. Observe that, by (|4.9[) we have 

V-C e (0)Wdxds = f f £^)|Vt?| 2 da;ds > e f f |W| 2 d:cds, 
io Jn Jo Jn Jo Jn 

and analogously for the term involving V£ e (i? s ). Therefore, we obtain 

1 f \L £ ($(t))\ 2 dx + e f ( \Vti\ 2 dxdr + \ ( \L £ (i) s {t))\ 2 dx + e f f |W s | 2 da;dr 

2 Jn Jo Jn 2 J Vo J Q J Vo 

f k{X){d-d s ){Z £ {d)~L £ {d s ))dxdr 
\ Jr * (4.98) 
t\§~V s )W E {uu)n)\\d t u T \{L e (V)-L e (V s ))dxdr 



o Jr c 

<^JjL £ ($ £ )\ 2 dx+±^ \L £ (^ S )\ 2 dx + I 12 + I 13 + h 4 . 



Now, due to the monotonicity of L £ the fourth term on the left-hand side is non-negative, and so is 
the fifth one, as 

c'(y - z)(L £ (y) - L s (z)) = c'(y - z)(y - g )*M^*M > for all y ^ z 

y-z 

also in view of (|3.9p . On the other hand, by the very definition (|4.2p of L £ , there holds 

\\L £ ($o\)\\ 2 h < ^olll +2|| M*o)lll < C (4.99) 
thanks to (|4. 14c|) - (|4.14d|) . and we have an analogous bound for ||£ e (^s' e )|| 2 j r • Moreover, we estimate 

h2 = - [ f div(d t u)L £ (tf)dxdr < f \\d t u\\ 2 v ds + ] f \\£ e (4)\\ 2 H dr 
Jo Jn Jo 4 Jo 

h 3 = - [ [ hL £ (#)dxdr< f \\h\\ H \\£e(4)\\Bdr 
Jo Jn Jo 

Im = - [ I \'(X)d t XL £ ($ s )dxds<C f (l + \\X\\ L oo ( r c ))\\d t X\\ Hr jL e (>& 8 )\\ Hrc dr. 



o Jr c 

We plug the above estimates into the r.h.s. of (|4.98[) . and use the previously proved bounds (|4.70l) . 
(|4.75p (yielding a bound for X in L 2 (0, T; L°°(L C )) ), and (|3.17l) for h. Applying a generalized version 
of the Gronwall Lemma (see, e.g., [6]), we conclude 

||£ e WIUo=(0,T;H) + ||£e(tf S )|U~ (0 ,T;tf rc ) < C ( 4 - 100 ) 

whence, in view of (14.681) . 

||ln e (^)|| iO o (0>x . ff) + ||ln e (tf s )|| L oo ( o, T; ff r j < C. (4.101) 
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Remark 4.12 As previously mentioned, the First, Fifth, Sixth estimates should be performed on a 
further approximate version of Problem (Pe). In fact, identities (|4.77|) . (|4.78j) . (|4.86|) . and (|4.94|) arc 
just formal since the 3t£ e (i?) only belongs to L 2 (0,T;V') (analogously, c?i£ £ (?9 s ) only belongs to 
L 2 (Q,T; (ff 1 (r c )) / ). These calculations can be rigorously justified in a framework where equations 
(|4.25p and (|4.26|) arc further regularized by adding viscosity contributions modulated by a second 
parameter v > , that is 



/ d t L e ($)vdx- [ div(9 t u)wda;+ / \Jfi\Jvdx + v [ V(d t d)Vvdx 
Jn Jn Jn Jn 



+ J k(X)(&-<& s )vdx + J c'(i?-i? s )|^(^(uN)n)ptUT|vda:= (h,v) v VveV a.e. in (0,T), 

(4.102) 

/ dt£ e (tf s )vdx - / d t X(X)vdx+ f S7ft s Vvdx + v f V(<9 t tf s ) Vvdx 

Jv c Jr c Jr c (4.103) 

k(X)(-8-ti s )vdx + c'(-&-# s )\$.(4>' e {uN)n)\\dtUT\vdx Vv£Vr c a.e. in(0,T), 



The presence of these additional viscosity terms in (|4.102p and (|4. 103[) implies that the so- 
lution to the PDE system of Problem (P e ), with (|4.102j) in place of (|4~2"5]l and (|4.103p in place of 
(I4.26p . and supplemented by natural initial conditions, satisfies in addition 

,C e (0) ) 0eff 1 (O,T;V), £ 8 (i?.),<?. 6fl" 1 (0,T;Vr e ) (4.104) 

and hence the formal identities (|4.77[) . (|4.78|) . (|4.86p . and (|4.94p can be rigorously revised. Such an 
approximation of Problem 13. 31 was considered in [3], (see also [TD] ) , to which we refer the reader. Let 
us just mention here that, for technical reasons (see [9j Remark 3.2]) the viscosity parameter v has 
to be kept distinct from Yosida parameter e for the logarithm. Hence it is necessary to derive global 
a priori estimates independent of e and/or v , and then perform the passage to the limit procedure 
in two steps, first as v \ and subsequently as e \. . 

To avoid overburdening the paper, we have preferred to omit this further vanishing viscosity 
regularization, at the price of developing the calculations for the First, Fifth, and Sixth estimates 
only on a formal level. 



5 Proof of Theorem [T] 

In this section, we detail the passage to the limit in the approximate Problem (P e ) as e tends to 
and we achieve the proof of Theorem [1] showing that the approximate solutions converge (up to 
a subsequence) to a solution of Problem 13.31 Hereafter we make explicit the dependence of the 
approximate solutions on the parameter e and use the place-holder 

Ve '■= ^K) n ' 

We split the proof in some steps. 



Compactness. Combining estimates l|4.68p ~ (|4.72p . (|4.75|> — (|4.T6Jt . and (|4.82|) with the Ascoli-Arzela 
theorem, the well-known (3UJ Thm. 4, Cor. 5], and standard weak and weak * -compactness results, 
we find that there exists an nine-uple w, 'd s , w s , vi, X, rj,z, £) such that, along a suitable (not 
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(5.1) 
(5.2) 



relabeled) subsequence, the following convergences hold 
u e ^u inH^OjTjW), 

u e ^u in C o ([0,T];j?' 1 -' 5 (n;]R 3 )) for aU 5 G (0, 1], 
X £ ^*X in L 2 (0,T;H 2 (T c )) n L°°(0,T;Vr e ) H -H" x (0, T; ifrj, 

X £ X in I 2 (0,T;ff 2 -' , (r c ))nC ([0,r];if 1 - i (r c )) for all p G (0,2] and S G (0,1], 

inL 2 (0,T ;j ffrJ, (5.3) 

T7 e -7j inL 2 (0,T;Y^), (5.4) 

z e ^*z in i°°(r c x (0,T);M 3 ), (5.5) 

jLt £ = |^(r; e )|z e ^*/ti in L°°(0,T;L 2+v {T c ;R 3 )) with > from (J3IT3J, (5.6) 

in L 2 (0,T;1/), ez? E -> in L°°(0, T; if), (5.7) 

tf s , e ^tf s ini 2 (0,T;FrJ, etf a , e -> in L°°(0, T; iT r J, (5.8) 
£ e (tf £ )^*u; in L°°(0, T; if) n if J (0, T; V') , 
£ e (tf £ )^™ inC°([0,T];F'), 
£e(#,,e)^*w a in i°° (0, T; ff r J nif^O, T;V/ e ), 
£ £ (tf s , £ )^™ s inC°([0,T];VS) 



(5.9) 
(5.10) 



as e J, . In addition, in view of condition (|3.13p on 31 , we have 

3?(rj e ) -> 3l(rj) in i°°(0, T; L 2+!/ (r c ; M 3 )), so that fj, = |#(tj)|z, (5.11) 
and (|5.6p improves to 

|K(»7 e )|z £ -* |K(»7)|z inL oo (0 > T;L a+, '(r e ;R 3 )). (5.12) 

Next, applying a generalized version of the Aubin-Lions theorem for the case of time derivatives as 
measures (see, e.g., [221 Chap. 7, Cor. 7.9]), from (|4.69[) . and estimates (|4.92[) and (|4.96[) we deduce 
that 



(5.13) 



tf £ ->tf ini 2 ^,^^ 1 -*^)) for all 5 G (0,1], 

ini 2 (0,T;L 5 (r c )) for all 5 G [1, 4), 
tf s , £ i? s in i 2 (0, T; L 5 (r c )) for all 8 G [1, +oo). (5.14) 



Moreover, taking into account the Lipschitz continuity and the C 1 - regularity of c (cf. (|3.9p ). from 
([5TT3]) - ([57Iill . we have 



c(z9 e -0 BtB ) -»- c(i?-0 a ) ini 2 (0,T;L 5 (r c )) for all tf G [1, 4) , 
c'(tf £ -tf s , £ ) ->V(0-0.) hiL q (0,T;L q (T c )) for all 9 G [1, oo) 



(5.15) 



Passage to the limit in (|3.44[) . Now, we consider (|3.44p - (|3.45[) written for the approximate 
solutions (u e , X e , $ Sj£ , £ e ) e . Taking into account convergences (|5.ip - (|5.3p . (|5.14p . and the Lipschitz 
continuity of A' and a' (cf. (|3.15l) . (|3.16I0 . we easily conclude that the limit quadruple (u,X, # s ,£) 
satishes equation (|3.44p . Combining the weak convergence (|5.3p with the strong one specified in (|5.2p , 
and taking into account the strong-weak closedness in L 2 (0,T; Hr c ) of the graph of (the operator 
induced by) (3, we conclude that £ £ (3(X) a.e. on T c x (0,T), i.e. (|3.45p holds. 

Passage to the limit in P~2"7]) . Owing to convergences ([53]) -- ([531) . (|5~4)l , (|5~6l) - (|5~7f and ([BTT5]) . 

we can pass to the limit in (|4.27|l . We get 

6(<9 t u,v)+a(u,v)+ / i?div(v)dx+ / Xuv dx+(r] ,v) Yrc + c(#-i? a )/z-vds = (F, v) w , (5.16) 
Ji2 7r n Jr, 
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for all v € W. Now we have to identify 77 and /j. as elements of d(p(u) and Tf)\d(dtu) , 
respectively, i.e. to show that (|3.42[) and (|3.43p hold. 

First, we test (|4.27|) by u e . For every t e [0, T] we have 

Um sup / / r] £ ■ u £ dx ds 

£-*o Jo Jr c 

= - liminf (b(u e (t), u E (i)) - 6(u , u ) + / (a(u e ,u e )+ / $ e div(u e ) dx) ds 
v Jo Jn 

+ [ (I X e \u e \ 2 dx+ [ c(i? £ - 1? S)E K • u £ dx) ds - / (F,u £ ), 
Jo Jr e Jr c Jo 

<- / ( b(d t u, u) + a(u, u) + / tfdiv(u)da: + / X|u| 2 dx+ / c(# -u- (F,u) w ) ds 

Jo V Jo Jr c Jr c 







T7 • u dx ds 



where the < follows from exploiting (|5.1|) . (|5.2p . (|5.4p . (|5.6|) . (|5.7p . and (|5.15l) . combined with lower 
scmicontinuity arguments, and the last equality is due to (|5.16|) . We use the above inequality and to 
show that for all v e Yr c and t € [0, T] there holds 



/ (77, v — u}Y rc ds < liminf / / T7 £ • (v — u e ) dxds < liminf / / 



j(un) - ^ £ (mn)) dxds 



< / / (0(i'n) - 0(un)) dsds 
Jo Jr c 

= / M v ) - ¥>(")) ds, 
Jo 

where the second inequality follows from the fact that rj e = (f>' £ (u^)n , the third one from the Mosco- 
convergence (see, e.g., [1]) of <f> e to <j>, and the last one from the definition f|3 . 1 2[) of ip . All in all, 
we conclude (|3.42p . 

Let us now show (j3~43)l . Preliminarily, for every fixed § G L 2 (0,T;V), d s G L 2 (0, T; Vr c ), 
and r] G L 2 (0,T;Y^ e ), we introduce the functional 3^ s ,rj) ■ £ 2 (0, T; L 4 (r c ; R 3 )) -> [0, +00) 
defined for all v G L 2 (0, T; L 4 (T C ; R 3 )) by 

3<fi,* m ,m (v) := / / c(i9(x,t)-^(x,t))|3?(7 7 )(a ; ,t)|j(v(x,t))dxdi 
Jo Jr c 

/ c(i?(a;,t)-^ s (ar ) t))|X(T7)(a;,t)||vT(a; ! t)|da:dt 
Jr c 

Clearly, 3($,'& 3 ,r)) is a convex and lower semicontinuous functional on L 2 (0, T; L 4 (r c ; R 3 )) . It can 
be easily verified that the subdifferential d3(-»,-» s ,rj) ■ L 2 (0 7 T; L 4 (r c ; R 3 )) =3 L 2 (0, T; L 4 / 3 (r c ; R 3 )) 
of ?(<>,<>„, rj) is given at every v G L 2 (0, T; L 4 (r c ; R 3 )) by 

fheL 2 (0,T;i 4 / 3 (r c ;R 3 )), 

heagMfl tjh v ) ^ ■( 5.17) 
1 \h(a: I t)Gc(<?(x > t)-^(a:,t))|31(T|)(ar,t)|d(v(a: I t)) 

for almost all (2, t) G r c x (0, T) , where d = 9j is given by (|1.13[) . We shall prove that 

3(*,t.,n)W -3 ( ^..rj)(9tu) > / / c(i9-i9 s )|K(tj)|z- (w-d t u)d:rdi (5.18) 

Jo Jr c 

for all w G L 2 (0, T; L 4 (r c ; R 3 )) . From (J535J we will conclude that c(i?-i? g )|3l(?7)|z € d2tf,#„n)(d t VL) , 
hence the desired (|3.43[) by (|5.17p . the strict positivity (|3.9p of c, and (|5.1ip . In order to show (|5.18l) . 
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we first observe that 



lim sup / / c(i? e -i? s , e )|^(77 e )|z £ -StUedxdt < / / c(i9-^ s )|3l(r;)|z-5 t udxdt, (5.19) 
e->o Jo Jr c Jo Jr c 

which can be checked by testing (|4.27p by d t u 6 and passing to the limit via convergences (|5.ip - (|5.2p . 
(|5.4p - (|5.6p . (|5.1ip . (|5.13p . lower semicontinuity arguments, and again the Mosco convergence of <j> e . 
Therefore, we have 

/ / c(#-i9 s )|DI(t/)|z- (w-d t u)dxdt <liminf / / c(i5» e - d s , E )\3l(r] e )\z e • (w - d t u e ) dxdt 
Jo Jr c Jo Jr c 

liminf/ / c^-^I^JKIwTl-K^TDdxdt 
Jo Jr c 



< 
< 



/ / c(i?-7? s )|^(77)|(|w T |-|(9tu) T |)dxdf (5.20) 
Jo Jr c 

where the first inequality follows from (|5 . 19[) and convergences (|5 . 1 2[) and f|5 . 15[) , the second one from 
the fact that \"R(ri £ )\z e <E \3l(r) £ )\d(dtu e ) , and the last one from combining the weak convergence 
(|5.ip with the strong convergences (|5.1ip and (|5.15p . Then, (|5 . 1 8[) ensues. Furthermore, arguing 
as in the derivation of (|5.20p . relying on (|5.ip . (|5.1ip . (|5.15p . and (|3.43p , and using that, indeed, 
\31(t] e )\z £ ■ d t u £ = ft(r/ e )||(<9 t u 6 ) T | a.e. in T c x (0, T), we deduce 

liminf / / c(i? e - i? s . e )|3?(r/ e )|z e • d t u e dxdt 
Jo Jr c 

= liminf / / c(d e -$ s . E )\X(r] e )\\(dtU e h\dxdt (5.21) 
Jo Jr c 

> f [ c(ti-# s )\$.(r])\\(d t u) T \dxdt= / / c(0 - 0,)|3l(»j)|ss • d t udxdt . 
Jo Jr c Jo Jr c 

Ultimately, from (|5.19[) and (|5.2ip we conclude 

lim / / c{d E -$ s ,)\jl(r] e )\z e -dtu E dxdt= / / c(tf - tf s )|rR(rj)|z • <9 { uda;dt. (5.22) 
<^°Jo Jr c Jo Jr c 

Now, in addition to (|5.ip . we prove the following strong convergence 

d t u £ ^d t u in L 2 (0,T;W), (5.23) 

which is crucial in order to pass to the limit in the frictional contribution / c'(# e — i? s 6 )|3?(</>t(u ejv n)| |(dtu e ):r|w 
in (|4.25p and (|4.26p . To this aim, we first observe that 

lim sup / b{d t u e7 d t u E )dt < / b(d t u,d t u)dt (5.24) 

e->0 Jo ^0 

arguing in a similar way as in the derivation of (|5 . 19|) : we test (|4.27p by dtu e and we pass to the 
limit exploiting convergences (|5.ip - (|5.2p , (|5.4p - (|5.6p . (|5.13p . (|5.22|) and the Mosco convergence of 4>e ■ 
Since the converse inequality for the lim inf e _>.q holds by the first of (|5.ip , we then conclude that 

lim / b(dtu £} dtU e ) dt = / b(dtU,dtu) dt. 
£ ~ >0 Jo Jo 

This gives |[5^5jl . by the W -ellipticity of 6 (cf. ([3T5])). 
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Passage to the limit in (|4.25[) and (|4.26j) . We pass to the limit in (|4.25p and (|4.26[) relying on the 
above convergences (f5Tj) - (f53|) . (fO|) . ff5T7 |) — (f5TTT] ) . (|5~T3l) - (|5~T5l) . (|533l) . and the following additional 
convergences for the nonlinear terms in (|4.25|) and (|4.26p . Indeed, conditions (j3.14p - Q3.15p on k and 
A and convergences f|5 . 2[) . (|5.13|) (|5.14l) yield 

k{X s )(# s - # s , £ ) ^ k(X)(ti - # s ) in L 2 (0, T; #r c ) , (5.25) 

as well as 

X(X e )^X(X) inF^T;!/ 5 / 2 ^)). (5.26) 
Exploiting all of the above convergences, we get 

(d t w,v) v — / div(<9 4 u) v dx + / V^Vudx 

Jn Jn (5.27) 

+ / k(X)('d-<& s )vdx+ / c'(i?-i? 5 )|D?(j7)||9 t u T |uda; = VueV a.e. in (0,T), 



(5.28) 



(<9 t w s ,z;) y - / d t X(X)vdx+ [ VdsVvdx 
Jr c Jr c 

= J k(X)(ti -tf s )vdx + J c'(tf -tf s )\R{r])\\dtu T \vdx \/v G Vr c a.e. in (0, T) . 



It remains to show that 

w(x,t) = ln(-d(x,t)) for a.a. (x,t) G x (0,T), (5.29) 

iy s (a:, i) = ln(tf s (a:, i)) for a.a. (x, t) G T c x (0, T) . (5.30) 

We argue just for (|5.29p . the procedure for (|5.30[) being completely analogous. First, recalling the 
definition of £ e (cf. (|4.2p ). we observe that (|5.9p and the second of (|5.7p give 

ln e (i? e )^*w inL°°(0,T;F). (5.31) 

Thus, by relying on well-known properties of Yosida regularizations (cf. Lemma 1.3, p. 42]), to 
conclude (|5.29|) it is sufficient to check that 



limsup [ I Ya. e {$ e )'d £ dxdt < [ ( wddxdt. 
e\o Jo Jn Jo Jn 



(5.32) 



The latter follows combining the weak convergence (|5.31[) with the strong convergence (|5.13|) . This 
concludes the proof. □ 
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